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ABSTRACT 

The Ovi ion observed in quasar absorption line spectra is the most accessible 
tracer of the cosmic metal distribution in the low redshift (z < 0.5) intergalactic 
medium (IGM) . We explore the nature and origin of O vi absorbers using cosmological 
hydrodynamic simulations including galactic outflows with a range of strengths. We 
consider the effects of ionization background variations, non-equilibrium ionization and 
cooling, uniform metallicity, and small-scale (sub-resolution) turbulence. Our main 
results are (1) IGM Ovi is predominantly photo-ionized with T w io 4 - 2±0 - 2 K. A 
key reason for this is that O vi absorbers preferentially trace over-enriched (by ~ x5) 
regions of the IGM at a given density, which enhances metal-line cooling such that 
absorbers can cool to photo-ionized temperatures within a Hubble time. As such, O vi 
is not a good tracer of the Warm-Hot Intergalactic Medium. (2) The predicted O vi 
properties fit observables if and only if sub-resolution turbulence is added, regardless 
of any other model variations. The required turbulence increases with Ovi absorber 
strength. Stronger absorbers arise from more recent outflows, so qualitatively this 
can be understood if IGM turbulence dissipates on the order of a Hubble time. The 
amount of turbulence is consistent with other examples of turbulence observed in the 
IGM and galactic halos. (3) Metals traced by Ovi and Hi do not trace exactly the 
same baryons, but reside in the same large-scale structure. Our simulations reproduce 
observed alignment statistics between Ovi and Hi, yet aligned absorbers typically 
have Ovi arising from cooler gas, and for stronger absorbers lower densities, than 
H i. Owing to peculiar velocities dominating the line structure, coincident absorption 
often arises from spatially distinct gas. (4) Photo-ionized O VI traces gas in a variety 
of environments, and is not directly associated with the nearest galaxy, though is 
typically nearest to ~ 0.1L* galaxies. Weaker Ovi components trace some of the 
oldest cosmic metals. (5) Very strong absorbers [EW > 100mA) are more likely to be 
collisionally ionized, tracing more recent enrichment (< 2 Gyr) within or near galactic 
halos. 

Key words: intergalactic medium, galaxies: evolution, cosmology: theory, methods: 
numerical 



1 INTRODUCTION 

The exploration of metals in the low redshift (z < 0.5) 
intergalactic medium (IGM) has taken a large leap for- 
ward with the advent of high-resolution space-based ul- 
traviolet (UV) spectroscopy with the Hubble Space Tele- 
scope Imaging Spectrograph (STIS) and Far Ultraviolet 
Spectroscopic Explore (FUSE). The strongest and most 
common metal transition seen in quasar absorption line 



cent papers (e.g 


. Tripp & Savagel 200rt Savage et al.l 20021: 


Prochaska et al. 


2004k iRichter et all 120041: ISembach et all 


2004al: [Lehner et al. 2006: Cookscv et al. 2008T) have exam- 



ined O vi absorbers along single sight lines, and have fit ion- 



ization models to individual systems with associated Lya 
and lower ionization metal transitions. Three recent stud- 
ies (Tripp et al. 2008, hereafter T08; Danforth & Shull 2008, 
hereafter DS08; and Thorn & Chen 2008a, 2008b) have com- 
piled samples of O vi, providing the largest statistical sur- 
veys to date. 

Understanding O vi in the low-z Universe is ex- 
tremely important, because it may hold the key to lo- 
cating a signifi cant r e servoi r of cosmic baryons and met- 
als (e.g. iTripp et al.l l2000lh The inventory of observed 
baryons (e.g. Fukugita et al.l Il998f ) falls well short of the 
predicted cos mological values from the cosmic microwave 
background (|Hinshaw et al.l l2008l h leading to the well- 
known "missing baryons problem." Cosmological Simula- 
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tions l|Cen fc Ostrikerl 1 19991 ; iDave et all Il999l . l200lh sug- 
gest that a large fraction of baryons today (> 30%) re- 
side in a hard-to-observe warm-hot intergalactic medium 
(WHIM), with temperatures of T = 10 5 - 10 7 K. The 
WHIM results primarily from shocks created during large- 
scale stru cture formation. More recent simulations includin 
feedback (|Cen fc Qst rikcr 20061; [oppenheimer fc Davefioo" 
hereafter OD08) find that galactic winds may increase the 
fraction of cosmic baryons in the WHIM to 50% or more. 

A possibly related problem is that o nly | of metals have 
been accounted for observationally (Fu kugita fc Peebles! 
120041 '). judging from the mismatch between the amount of 
metals nucleosynthesized and ejected by observed stars and 
those observed in various cosmic baryonic components. O VI 
potentially represents the Holy Grail of all things missing 
in the low-z Universe, because its collisionally ionized equi- 
librium (CIE) maximum temperature, 10 5 ' 45 K, provides a 
unique and relatively easily accessible tracer of WHIM met- 
als and baryons. The incidence of Ovi with broad Lya ab- 
sorbers (BLAs) th ought to trace gas at T > 10° K (e.g. 
iRichter et al.l l2004h supports this notion, while O VI seen 
with very broad Hi having line widths > 1 00 kms~ may 
indic ate even hotter O VI at T ~ 10 6 K (e.g. iDanforth et al.l 
120061 ). 

Earl y investigations i nto Ovi using cosmological sim- 
ulatio ns (ICen et al l I200II ; iFang fc Brvanll200ll ; IChen et all 
120031 ) predicted that stronger O VI absorbers tend to be col- 
lisionally ionized while weaker ones tend to be photo-ionized, 
with the cross-ove r equivalent width of ~ 30 — 50 mA. 
ICen fc Fan e (2006) added non-equilibrium ionization and 
used higher resolution simulation, finding similar behavior 
compared to their earlier work. The above simulations gener- 
ally focused on fitting only the observed cumulative equiva- 
lent width (EW) distribution. Recent observations now pro- 
vide new challenges for simulations to fit a wider range of 
low- z Ovi observables. 

Recent surveys have renewed confusion about the na- 
ture of O VI absorbers. Ionization models for absorbers show- 
ing aligned Ovi, Hi, and Cm are usually forced to invoke 
multi-phase gas, with both O VI in CIE a nd lower ioniza- 
tion species at photo- ionized temperatures (IProchaska et al.l 
12004 IDanforth etHI 120061 ; ICooksev et al.ll2008l ). Yet T08 
find that the majority of aligned H I — O VI absorbers can 
also be explained by photo-ion ization alone. The "bl ind" 
searches of O VI without H I by iThom fc Cher] (1200881 ) in- 
dicate > 95% of O VI is associated with H I. These stud- 
ies suggest that Ovi traces baryons at least partially ac- 
counted for by H I absorbers, and is not necessarily a good 
tracer of the WHIM. Studies of galaxy-absorber correla- 
tions find a wide variety of environments for Ovi ab- 
sorbers including v oids, filaments, galaxies , and groups (e.g. 
Tripp et all 1200 it ; iTumlinson et all 120051 ; IProchaska et all 



20061 ; iTripp et all I2006T T Furthermore, it is not clear how 



the Ovi in intermediate- and high- velocity clouds (IVCs 
and HVCs) associated with the Milky Way (MW) halo (e.g. 
ISembach et~alll2003l ; ISavage et aill2003l ; iFox et al.ll2005h re- 
late to O VI observed in quas ar absorption line spe ctra. A 
unified model was proposed bv lHeckman et all (|2002T ) where 
O VI absorbers are all radiative cooling flows passing through 
the coronal temperatu re regime, with a modification by 
iFurlanetto et al. (2005) for long cooling times in the IGM. 
However, observations are inconclusive as to whether this 



scenario also appli es to IGM Ovi (e.g IDanforth et al .1120061 : 
iLehner et all 120061 ). Hence there remains much uncertainty 
in regards to which cosmic gas and metal phases Ovi ab- 
sorbers actually trace. 

In this paper we explore Ovi and Hi absorbers in a 
range of physical models using our modified version of the 
cosmological hydrodynamics code Gadget-2. We pay close 
attention to three key observables confirmed in multiple 
studies: the cumulative EW distribution, O VI linewidths (6- 
parameters) as a function of O VI column density, and the 
alignment of Ovi with Hi. Our two main goals are (1) to 
see how self-consistent metal enrichment via galactic super- 
winds reproduces the O VI observations, including an honest 
evaluation of the short-comings of our simulations, and (2) 
to understand the physical conditions and environments of 
O VI absorbers. This study combines state-of-the-art mod- 
eling with the most recent O VI data to further understand 
the nature of IGM O VI, in light of the much anticipated 
installation of the Cosmic Origins Spectrograph (COS) and 
the re-activation of STIS on Hubble. 

The structure of the paper flows as follows. In §2, we in- 
troduce our simulations all run to z — with three different 
galactic outflow models, plus other post-run input physics 
variations. Then, in §3 we see how the various models fit the 
Lya forest and our three O VI observables, including a test 
of resolution convergence. The crux of this paper lies in the 
following sections where we dissect our simulated absorber 
population, acknowledging the imperfections in our model- 
ing while assessing what we believe is physically significant. 
§4 advocates that IGM O VI absorbers are primarily photo- 
ionized, with an in-depth analysis of the association with H I, 
and forwards an explanation for the non-thermal component 
of Ovi line profiles as arising from small-scale turbulence. 
§5 examines the origin of Ovi absorbers with an emphasis 
on environment. We discuss the minority population of colli- 
sionally ionized O VI here, and examine tell-tale signs of such 
WHIM absorbers with an eye towards COS. We summarize 
in §6. 

Throughout this paper, we adopt [Xsplund et all (|2005l ) 
solar abundances, as most low- z O VI observations use 
this or a similar value. This contrasts to our use of 
lAnders fc Grevessel (|l989l ) values in previous publications 
(Oppenheimer & Dave 2006, hereafter OD06; Dave & Op- 
penheimer 2007; OD08) , and results in a significant decrease 
in solar oxygen mass fraction (0.00541 versus 0.00962). The 
Type II s upernovae (SNe) yields w e use as simulation inputs 
are from iLimongi fc Chieffil l|2005h . and are independent of 
solar values. For ease of discussion, we classify O VI absorbers 
into three categories according to their column density; 
weak (TV(Ovi) < 10 13 ' 5 cm" 2 ), intermediate (TV(Ovi) = 



10 1 



cm ), and strong (iV(0 vi) 



10 1 



We often split the intermediate absorbers into two bins, since 
this comprises the majority of observed O VI. 



2 SIMULATIONS 

We use cosmological simulations run to 2 = using our mod- 
ified version of the N-bo dy + Smoothe d Particle Hydrody- 
namics code Gadget-2 (Springel 2005) in order to explore 
the nature of Ovi in the low- 2 Universe. Our simulations 
directly account for cosmic metal enrichment via enriched 
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galactic outflows, where outflow parameters are tied to 
galaxy properties following that observed in local starbursts. 
Our most successful model follows scalings for momentum- 
driven winds; this fairly uniquely matches a wide range of 
observations such as early IGM enrichment (OD06), th e 
galaxy mass-metallicity relations (|Finlator fc Dave] [2008), 
and the enrichment a nd entropy levels seen in intragroup 
gas (|Dave et al.ll2008T h We introduce the simulations first, 
then describe variations applied to our favored outflow 
model, and finally describe post-run variations to the input 
physics. 



2.1 Model Runs 

Our simulations adopt cosmological parame ters based on 
5-year WMAP results (|ffinshaw et all 120081 ). The param- 
eters are Q = 0.25, Q A = 0.75, f2 6 = 0.044, H = 70 
km s _1 Mpc -1 , <T8 = 0.83, and n = 0.95; we refer to 
this as the d-series. The value of as is slightly higher 
than favored from WMAP data alone, but it agrees bet- 
ter with combined results from Type la SNe and baryon 
acoustic oscillation data. Our general naming convention is 
A[boxsize]n[particles/side] [wind mode!]. 

We use the momentu m-driven wind model based on 
the a nalytic derivation by iMurrav. Quatert. fc Thompson! 
l|2005l) and explored previously in OD06 and OD08. Briefly, 
wind velocity (ivtnd) scales linearly with the galaxy velocity 
dispersion (a), and mass loading factor (77) scales inversely 
with a. The mass loading factor represents the mass loss 
rate in outflows in units of the galaxy star formation rate. 
We use the following relations for u w ind and 77: 



800 



Wwind = 3<T y/zT 
(To 

n = — , 



1, 



(i) 

(2) 



where fi is the luminosity factor in units of the galactic Ed- 
dington luminosity (i.e. the critical luminosity necessary to 
expel gas from the galaxy potential), and o~o is the normal- 
ization of the mass loading factor. A minor difference from 
OD06 and OD08 is that we no longer impose an upper limit 
on v w ind due to SN energy limitations. 

The model runs, detailed in Table [T] have [wind model] 
suffixes denoting a variation of the momentum-driven wind 
model run. Our previously used favored model, vzw, contains 
a spread in fi — 1.05 — 2.00, includes a metallicity depen- 
dence for fi owing to more UV photons output by lower- 
metallicity stellar populations, and adds an extra boost to 
get out of the galaxy potential well simulating continuous 
acceleration by UV photons. 

Figure [1] illustrates a 25 kms -1 slice spanning the full 
extent of the d32n256vzwl50 simulation box at z = 0.25. 
The location of galaxies (upper left) is well correlated with 
the enriched IGM (upper right). The Lya forest traces 
baryons in filaments (bottom left) down to the lowest observ- 
able column densities. The Ovi absorbers (bottom right) 
mainly trace the diffuse IGM enriched via feedback, show- 
ing a more grainy structure indicative of inhomogeneous en- 
richment. Much of this paper focuses on quantifying the re- 
lationship between these various observational and physical 
components of the IGM. 

We run 16 and 32 /i _1 Mpc (comoving) sized boxes at 
the same mass and spatial resolution (with different numbers 
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Figure 2. Wind properties as a function of redshift for 
our three wind models, plus one higher resolution simulation 
(dl6n256vzwl50), and a simulation to test box size convergence 
(dl6nl28vzwl50). In order to create a reasonable baryonic frac- 
tion in stars at z = 0, weaker winds (top panel) must have a 
higher mass loading factor (middle panel) resulting in greater en- 
ergy input per unit star formation as shown by the E wind / Eg^ 
ratio (lower panel). The opposite applies for the Izw model. E$n 
is set to 1.3 X 10 49 ergs M" 1 . 



of particles) of this wind model. The baryon fraction in stars 
is identical, while the wind speed and mass loading shown 
in Figure [2] are in good agreements from z = 6 — > 0, despite 
a volume difference of 8x (compare green vs. cyan lines). 
E w i n d/EsN decreases slightly in the smaller volume, because 
fewer massive galaxies having higher wind energy efficiencies 
form at early times. 

We also run a 16/i _1 Mpc box length, 2 x 256 3 parti- 
cle simulation to z — (dl6n256vzwl50) to explore resolu- 
tion convergence. This run with its smaller time stepping is 
the most computationally expensive run, using over 30,000 
CPU hours on Harpertown Intel Xeon processors. The wind 
properties in Figure [2] indicate lower ^wind's and higher rj's 
due to increased contributions from small galaxies below 
the resolution limit of the d32n256vzwl50 box. The stel- 
lar baryon fraction is 4% higher in this simulations, mainly 
due to late time star formation. We suspect this is due to 
the non-cosmological scale of this box, which cannot form as 
many groups and clusters at late times where star formation 
is inhibited, as well as the fact that it probes further down 
the relatively steep mass function of star-forming galaxies. 

In this paper we explore two other variations of the 
momentum-driven wind model that bookend, at high and 
low wind speeds, physically plausible winds that are able to 
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Table 1. Simulations 



Name" L b e c mg PH Wind Model erg n*/fij,(* = 0)^ 



d32n256vzwl50 


32 


2.5 


34.6 


Momentum-driven 


150 


0.102 


dl6nl28vzwl50 


16 


1.5 


34.6 


u 


150 


0.144 


dl6n256vzwl50 


16 


1.25 


4.31 




150 


0.144 


d32n256hzw075 


32 


2.5 


34.6 


Hybrid Momentum-SNe 


75 


0.061 


d32n2561zw400 


32 


2.5 


34.6 


Low-i> w ; n( j Momentum 


400 


0.116 



a vzw suffix refers to the momentum-driven winds, hzw suffix refers to momentum-driven winds with an added boost 
from SNe. Izw suffix refers to momentum-driven winds launched with 51% as much kinetic energy relative to vzw. 
All models are run to z = 

b Box length of cubic volume, in comoving h~ J Mpc. 

c Equivalent Plummer gravitational softening length, in comoving h~ 1 kpc. 
''Masses quoted in units of 10 6 Mq. 

e Normalization for the momentum-driven wind mass-loading factor where r] = co/cr, units in kms -1 . 
Fraction of baryons in stars. 




Figure 1. A 25 kms" 1 slice of IGM spanning the 32 comoving h 1 Mpcs of the d32n256vzwl50 simulation at z = 0.25. Upper left: The 
locations of galaxies are shown as colored points, where the color corresponds to stellar mass, overlaid on the gas overdensity represented 
in greyscale. Upper right: Color indicates the enrichment level; 90% of the IGM remains unenriched as the superwi nds are unable to 
enrich voids. Bottom panels: Hi (left) and Ovi (right) column densities in this 25 kms -1 slice, assuming a uniform lHaardt h Madaul 
(2001) ionization background. Ovi is not as extended as Hi, partially due to the small volume filling factor of metals, however the 
strongest Ovi absorbers are often found associated with sub-M* galaxies, unlike Hi, which has its strongest absorption around the large 
group in the center. The vast majority of O VI in our simulations is photo-ionized tracing overdensities of 10-200. 



enrich the IGM as observed. In addition to varying iVind, 
we must also vary oo to matc h the observed b aryon fraction 
in stars today of ~ 6 - 10% (|Cole et al.ll200ll) . 

The high-v w i n d case is denoted hzw in which a portion of 
the energy from Type II SNe is used to increase w w ind, while 
keeping the momentum-driven wind derived r\ the same. We 



use 1.9 x 10 48 erg Mg 1 , or one-seventh of the total supernova 
energy from a Chabrier IMF, assuming each star between 8- 
100 M Q explodes with 10 51 ergs. Additionally, we add 1.9 x 
10 48 erg Mg 1 of SNe energy to the thermal energy of each 
wind particle, and allow adiabatic expansion (i.e. no cooling) 
while the wind particle is hydrodynamically decoupled. 
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The hzw model is designed to explore the case where 
we push feedback to extreme physical limits with the spe- 
cific purpose to see if such strong feedback can make IGM 
metals both more widely distributed and hotter in the low- 
z Universe. In order to obtain a reasonable stellar density 
by \ow-z we set ao = 75 kms -1 , which results in 6.1% of 
baryons in stars at z — 0. Figure [2] shows the wind param- 
eters as a function of redshift. Ironically, the lower r\ of the 
hzw model reduces the efficiency of winds as indicated by the 
lowest E w i n d/ Esn ratio of all three wind models explored. 

To explore the low-« w i n d case, we run the Izw model in 
which the wind velocity is reduced by 50% relative to the vzw 
model. In order to prevent the formation of too many stars 
by z — 0, ao is set to 400 kms" , resulting in f2*/f2(, = 0.103 
at z = 0. The Izw model generates the most wind energy 
per baryonic mass {E w ind/ Esn), because of the high mass 
loading factor. 

2.2 Post-Run Physics Variations 

A variety of post-run treatments are performed on the model 
runs, using modifications within our quasar absorption line 
spectral generator, specexbin. Some of these variations are 
used to explore the consequences of unlikely but intuitive 
scenarios (i.e. a uniform metallicity distribution), while oth- 
ers are meant to account for physics not included in our 
simulations (i.e. non-equilibrium ionization). We list these 
variations below, and denote such model variations through- 
out the text with an additional suffix at the end of the model 
name. 

Z01: This model variation uses a uniform metallicity of 

0. 1 Zq for every SPH particle instead of the metallicity from 
the simulation. This unrealistic model presents an exercise 
to see how O VI observations vary if metals are uniformly dis- 
tributed. 0.1 Zq is the fiducial value that was needed in early 
simulation studies to match the equivalent width distribu- 
tion of low-z Q VI abso rbers l|Cen et al.ll200ll ; lFang fc Brvanl 
l200ll ; IChen etUl2003l ). It is also similar to the average z = 
oxygen content of the d32n256vzwl50 simulation, 0.13 Zq. 
We do not alter the gas densities or temperatures, so the 
ionization conditions remain the same. 

ibkgd: By defau lt, we use a spatially uni form ionization 
background given of lHaardt fc Madau (2001) background, 
with a contribution from quasars and star forming (SF) 
galaxies (10% escape fractio n) included. We use the n ormal- 
ization taken directly from lHaardt fc Madaul (|200ll ) with 
no adjustment, since as we will show t his reproduces z ~ 
Lya forest statistics quite well (see also lPaschos et aT]|2008t ) . 
The hydrogen photo-ionization rate is T — 1.84 x 10 _13 s _1 
at z = 0.25. To explore variations in photo-ionization, we 
consider two other uniform ionizat ion backgrounds: 

1. ibkgd-Ql uses intensity from the lHaardt fc Madaul (|200ll ) 
background, minus the contribution from SF galaxies. This 
mainly affects H I properties as SF galaxies contribute ~ | of 
the ionization at the Lyman limit, therefore the absorption 
from the Lya forest will grow. The ionization intensity at 
the Ovi ionization potential, 8.4 Rydbergs, remains mostly 
unaffected by subtracting star forming galaxies. This back- 
ground may be applicable in regions where the mean free 
paths of H I ionizing photons are shorter than O VI ionizing 
photons, e.g. denser regions in and aroun d galactic haloes. 

2. ibkgd- Q3 uses triple the quasar-only lHaardt fc Madaul 



|200ll ) background, which better reproduces the global 
statistics of the Lya forest. There are now more pho- 
tons at 8.4 Rydberg, which leads to photo-ionized O VI at 
higher overdensities. We explore this alternative since lit- 
tle is known about the true ionization field at these high 
energies. 

noneq: Our default case assumes ionization equi- 
librium for all species , wh ich we calculate using 
CLOUDY i|Ferland et al.l Il998h . To test sensitivity to 
this assumption for collisionally ionize d Ovi, we apply 
the no n-equilibrium ionization tables of lGnat fc Sternberg! 
(2007) as a function of temperature and metallicity 
rather than our C LOUD Y-generated equilibrium tables. 
iGnat fc Sternberg) (|2007f ) calculated non-equilibrium ion- 
ization fractions as a function of metallicity assuming 
collisional ionization only. High-metallicity shocked gas 
cools more rapidly than it can recombine, resulting in higher 
ionization states than expected at lower temperatures when 
assuming CIE. This has implications for O VI at densities 
where it is collisionally ionized, because ionization fractions 
of a few percent for Ovi can exist at well below its CIE 
peak at 10 5 ' 45 K if the gas is Zq or above. There does not 
exist a non-equilibrium treatment for photo-ionized gas as 
of yet, since this depends additionally on density; however 
these ionization fractions are more likely to be regulated by 
the intensity of the ionization background. Ideally it would 
be best to directly i nclude non-equilibrium ionization as 
m I (Jen fc Fand {2006), but we leave that for future work, 
and consider this test as bracketing the plausible impact of 
non-equilibrium ionization. 

To more fully explore the impact of non-equilibrium 
ion ization, we also run a d !6nl28hzw075 simulation with 
the lGnat fc Sternberg! l|2007l ) cooling tables as a function of 
metallicity. The metal line cooling peaks for individual ions 
are spread out over larger temperature ranges as individual 
ion coolants exist over a broader swath of temperature. The 
cooling is also of the order 2-4 x less in some places than 
ISutherland fc Dopital (Il993l ). because ion coolants exist at 
lower temperatures where they collide less rapidly. This run 
is denoted as dl6nl28hzw075-noneq. To ensure that this vol- 
ume is not too small, we run a dl6nl28vzwl50 simulation 
and find general convergence with d32n256vzwl50 for the 
star forming, wind, and IGM enrichment properties; there- 
fore this box size should be sufficient for comparison with 
larger volumes. 

bturb: The most critical post-run modification we em- 
ploy is the by-hand addition of sub-resolution turbulence. 
As we will show, our simulations do not reproduce the 
broadest observed absorption lines, especially those with 
iV(Ovi) > 10 14 cm -2 . Often, where observations see one 
broad component, our simulations would show a number 
of narrow components in a system. This is understand- 
able when one considers our simulations have a mass res- 
olution of mspH ~ 3.5 x 10 7 A/q, and meta ls may exist 
in structures of significantly smaller scale (e.glRauch et alj 
19991; ISimcoe et al.ll2006l ; iFrank et al.ll2007l ; lHao et alj|2007l 



Schave et al.ll2007l ). Even at the best-available instrumental 



resolution of STIS, thermal broadening at CIE O VI temper- 
atures cannot explain the bulk of the line widths in the T08 
dataset. Furthermore, the alignment with H I suggests that 
some non-thermal broadening mechanism is needed. Tur- 
bulence is often used as the overarching explanation (e.g. 
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iThom fc Chenll2008bl ). We will argue that Ovi absorbers, 
especially strong ones, are likely to be made up of numerous 
smaller metal concentrations with a range of velocities re- 
sulting in a component dominated by turbulent broadening. 

We add turbulence residing on scales below the SPH 
particle resolution ("sub-SPH turbulence") by adding a tur- 
bulent 6-parameter as a function of hydrogen density (uh)- 
With no physical guide available, we add just enough tur- 
bulence to approximately match the observed line width 
distribution as a function of column density. Specifically, 
we apply a linear fit to the STIS-only sample of T08 
fe(Ovi) - JV(Ovi) relation between 10 12 ' 9 - 10 14 9 cm' 2 
(b oba = 21.40 x log[AT(Ovi)] - 268.4 kms" 1 ), and subtract 
our relation from the d32n256vzwl50 lines of sight at STIS 
resolution (b no turb = 1-50 x log[iV(0 vi)] —8.9 kms -1 ), using 
the equation 

^turb ^obs ^rtoturb- (3) 

bturb is the turbulent broadening describing the small-scale 
motions within an SPH particle for this specific simulation 
resolution. We cast b tur b in terms of hh by employing the 
relation 

log[ra H ] = 0.66 x log [AT (O vi)] - 13.98 cm" 3 , (4) 

which is a dependence we find in i)4.1l Putting it together, 
the relation we use in specexbin is 

b tUTb = ^/l4051og[n H ] 2 + 156741og[n H ] + 43610 kms" 1 , (5) 

applied only over the range hh = 10 -5 ' 31 — 10 -4 ' 5 cm -3 , 
describing the majority of O vi within filaments but outside 
galactic halos. 

Since the sub-SPH turbulence from Equation [S] would 
extrapolate to artificially high values at higher densities, 
we use two more observational constraints from the O vi 
surveys of MW IVCs and HVCs to des cribe sub-SPH 
turbu lence associated with galactic halos. ISembach et al.l 
|2003l ) finds high- velocity O vi features have 6(0 vi) = 
40 ± 14 kms -1 , and suggest this O vi arises between the 
interfaces of cool/warm clouds and the T > 10 6 K galac- 
tic corona/intragroup medium at R > 70 kpc with n < 
10 -4 — 10 -5 cm -3 . We assume that these lines are pre- 
dominantly turbulently broadened, as temperature and spa- 
tial explanations fall well short, requiring another broaden- 
ing mechanism possib ly related to the motion of the gas 
|Sembach et al.l 120031 ). If the absorption arises from CIE 
Ovi (T ~ 10 5 - 5 K) then n H « 10" 45 cm -3 applies for 
this gas assuming approximate pressure equilibrium with the 
galactic corona. Therefore b tU rb ~ 40 kms -1 corresponds to 
nn = 10~ 4 5 cm" 2 both in Equation \5\ and in HVCs. 

ISavage et ail (|2003l ) find 6(Ovi) « 60 kms -1 in IVCs, 
which they find most likely reside in the galactic thick disk 
with O vi densities tracing tih ~ 10 -3 kms -1 for O vi flows. 
Using a linear interpolation with these two constraints, we 
find the relation 

bturb = 13.93 log(n H ) + 101.8 kms -1 (6) 

for 7ih = 10 -4 ' 5 — 10 _3 '° cm -3 , and assume a maxi- 
mum bturb = 60 kms" 1 for higher densities. Very rarely 
do our random sight lines intersect densities as high as 
nu = 10 -3 cm -2 , but such dense regions can create a strong 
absorption profile. 

At z = 0.25, the sub-SPH turbulent broadening is 



bturb = 13,22,40, and 51 kms" 1 at p/p = 20,32,100, and 
320 respectively. Besides producing wider lines, it turns out 
to also help explain various other observed Ovi absorber 
properties, as we show in <14.4I 



3 SIMULATED VS. OBSERVED Ovi AND Hi 
ABSORBERS 

In this section we place constraints on our models from ob- 
servations of the low-z Lya and O vi forests. We begin with 
the Lya forest, considering column density and line width 
distributions, and the evolution from z = 1.5 — ► 0. We 
then focus on three key Ovi observables: (1) the cumula- 
tive equivalent width distribution, (2) fo(Ovi) as a function 
of AT(Ovi), and (3) the AT(Ovi) as a function of AT (Hi). We 
explore how each of our models fits these observables, and 
explain why models the turbulent broadening provide the 
best fit. We end the section discussing numerical resolution 
convergence among our simulations. 

We generate 30 continuous lines of sight from z = 1 — > 
shot at a variety of angles through each simulation using 
specexbin. We use a continuum-normalized spectral tem- 
plate with 0.02 A bins convolved with the STIS instrumen- 
tal resolution (7 kms - ), and add Gaussian noise with 
S/N — 10 per pixel. This creates a large sample of high 
quality spectra covering Az = 15, larger than observational 
samples, with the intention to extract subtle differences be- 
tween our various models. For the sake of simplicity, we gen- 
erate H i-only and O vi-only spectra using only the strongest 
transition (i.e. Ly-a and the 1032 A line). See §2.5 of OD06 
for a mo re detailed descrip tion of specexbin. The AutoVP 
package |Dave et al.l ll997T ) is used to fit all component ab- 
sorption lines, yielding e.g. a sample of 475 O vi components 
with A r (Ovi) ^ 10 13 cm" 2 in the d32n256vzwl50-bturb 
model. 

3.1 H I observables 

The observationally well-characterized Lya forest provides 
baseline constraints on our simulations. A successful model 
of O vi absorbers must first reproduce observations of the 
low- z Lya forest, especiall y given that Ovi is us ually asso- 
ciated with Hi absorption (|Thom fc Chenll2008al ). The Lya 
forest traces the warm photo-ionized IGM ( i.e. T ~ 10 4 
K), w hich contains ~ 30% of cosmic baryons l|Penton et al.l 
l2004h and likely traces some of the same gas as photo-ionized 
O VI absorbers. We examine briefly three key H I observables 
spanning a range of physical and evolution properties that 
our simulations must match. 

Figure |3] shows the column density distribution (left), 
evolution of line densities from z = 1 — * (middle), and 
the N — b relation (right) from observations and four of 
our simulations, namely d32n256vzwl50, dl6n256vzwl50, 
d32n256vzwl50-bturb, d32n256hzw075. We display only 
these models because there exist only minor differences in 
the Lya forest amongst the other models. The comparison 
to data shows reasonably good agreement, except for an un- 
derestimate by at least a factor of two of strong Hi lines 
(EW > 240 mA). There are two likely reasons for this: The 
first is that the observations compared to have lower signal- 
to- noise and/or resolution than our simulated spectra, so 
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line blending may artificially enhance the observed equiva- 
lent widths; and second, many of these lines arise from gas 
within galactic haloes, where we may not be resolving halo 
substructures adequately. We note however that the higher 
resolution in the dl6n256vzwl50 simulation does little to 
resolve this discrepancy, so either we need much higher res- 
olution or else line blending issues are important. 

We find the low-z Lya forest is relatively unaf- 
fected by wind strength ; this is similar to the findings of 
iBertone fe White! (2006) at high-z, who find galactic winds 
do very little to alter the main statistical properties of the 
z ~ 3 Lya forest. The only difference worth noting is that 
the strong winds of the hzw reduce the frequency of strong 
Hi lines (EW > 240 mA) by about ~20%. 

We now compare our baseline vzw model with and 
without turbulence, since we will later show that the 
d32n256vzwl50-bturb model best fits the O VI observables, 
and we want to quantify how sub-SPH turbulence influ- 
ences the Lya forest. Ideally, turbulence should not affect 
Lya observables significantly since the low-z Lya-forest is 
well-d escribed by simul ations without any added turbu- 
lence (|Dave et all Il999h . and the 6-parameters are well- 
described primarily by thermal and Hubble broadening 
l|Dave fc Tripdl200ll ). 

The left panel shows the differential H I column density 
distribution, i.e. the number of absorbers per Hi column 
density interval per redshift pathlength d 2 n/dN(H l)dX, 



Penton et al 



from o ur models compared to power law fits by 
20041 . dashed line, z = 0.002 - 0.069) and 
2007, dotted line, z < 0.4). The different models show lit- 



Lehner et al 



tie variation from each other below iV(Hi) = 10 14 cm 2 , 
agreeing well with the slightly steeper power law of the 
ILehner et all (|2007T ) fit (d 2 n/dN(Rl) oc A^Hl)" 1 ' 83 ) over 
the same redshift interval. We use fits to Ly/3 for Af(Hl) ^ 
10 14 ' 4 cm -2 absorbers, but find no statistical difference ver- 
sus using only Lya for the frequency of strong lines fit with 
AutoVP. We find fewer higher column density lines than ob- 
served (as also indicated by our comparison to the dN/dz of 
high equivalent width lines), but we note that an observer 
typically climbs up the Lyman series to fit the strongest lines 
resulting in a more accurate fit. 

We consider the redshift evolution of Lya absorbers in 
the central panel of Figure [3] The evolution is rapid down 
to z ~ 2 but slows significantly below z ~ 1 (iBahcall et al.l 
Il99ll ; llmpev et all 1 19961 ; IWevmann et al l Il998l ). due to a 
decreasing ionization field strength counter-balancing de- 
creasing rates of recombination due to Hubble expansion 
|Dave et al.lll999T l. Matching this trend is important, be- 
cause the lack of evolution in Hi may well be reflected in 
O VI. In this paper we do not explore O VI evolutionary 
trends, leaving it for future work, but we broadly find a 
comparable amount of evolution as in H I. We count the fre- 
quency of lines, dN/dz, for intermediate (240 > EW ^ 60 
mA, AT(Hl) = 10 131-140 cm" 2 ) and strong absorbers 
(EW 240 mA, iV(Hl) = lo 14 " 17 cm" 2 ). Our simulated 
in termediate absorbers a gree with the evolutio n from higher 
z llJanknecht et al.l2006l ) to the local Universe (|Penton et al.l 
12004 ) . We do not calibrate the factor of two mismatch in 
stronger absorbers, because this more likely has to do with 
unresolved halo substructure and the breakdown of the op- 
tically thin assumption of the ionization background in such 



regions; this is beyond the scope of this paper, and is not 
very relevant to the bulk of O VI absorbers, as it turns out. 

The rightmost panel of Figure [3] shows the 6(H l)-iV(H i) 
correlation, with 1-a d ispersions, compiled from a number 
of quasar sight lines bv lLehner et alj (|2007| ). Our absorbers 
show a slow but steady increase until iV(Hl) = 10 14 cm -2 . 
Above iV(H i) > 10 14 ' 4 cm -2 , fitting to Ly/3 rather than Lya 
avoids saturation effects, and yields systematically smaller 
6-parameters in AutoVP. The ILehner et al] (120071) sample 
combines FUSE and STIS data, and could result in differ- 
ent 6-parameters than our uniform simulated STIS spectra 
(with typically higher S/N). The turbulent broadening starts 
making an impact in Hi lines with AT(Hi) > 10 140 cm -2 , 
tracing overdensities around 30, although it is not obvious 
whether the agreement with data is improved or worsened. 
A detailed comparison must wait for when the simulated 
spectra are generated and analyzed identically to observa- 
tions. 

An important marker of WHIM gas is broad Lya ab- 
sorbers (BLAs). Ideally, we would like to derive the fre- 
quency of BLAs since they are often thought to trace 
the same WHIM gas as collisionally ionized O VI (e.g. 
iRichter et al.l 120041 ) . However, continuum fitting uncertain- 
ties, spectral resolution, and S/N sensitivity likely make any 
such comparison fraught with systematics. We find BLA 
frequencies of dn(BLA) / dz — 16 and 12 using S/N = 20 
and 10, respectively, for AT (Hi) > 10 13 ' 2 cm -2 in the 
d32n256vzw!50-btu rb simulation, compared to 30 ± 4 from 
ILehner et~ al. (2007). However, 39% of these absorbers have 
iV(Hl) > 10 14 cm -2 , while the data shows a lower fraction 
of strong BLAs. We appear to be missing a population of 
weak BLAs; however we leave a detailed comparison of such 
properties for future work during the COS era. 

Our review of H I finds we reproduce most key observ- 
ables; only subtle disagreements among our various mod- 
els are noticeable, and then only for strong lines above 
AT (Hi) = 10 14 cm -2 . These disagreements may be related 
to (1) detailed differences between processing observed and 
simulated spectra, and (2) the inability of simulations to 
resolve absorbing substructures in halo gas. The strongest 
winds make 20% fewer AT (Hi) = io 14 °- 170 cm" 2 absorbers 
at all redshifts below z = 1, while added turbulence appre- 
ciably broadens only absorbers with A r (Hl) 10 14 cm -2 . 
Much remains to be learned about the observed low-z Lya 
forest by comparing simulations to future observations with 
COS, especially for BLAs; we leave for future work exam- 
ining the entire Lyman series with identical reduction and 
analysis methods applied to observed and simulated spectra. 



3.2 O VI Observables 

We now consider three robust O VI observables with trends 
corroborated by multiple groups. We explain their physical 
significance, and then explore how each of our models, all 
at the same resolution, fits these three sets of data. Lastly, 
we compare models at different resolutions to examine the 
issue of numerical resolution convergence. 

The cumulative EW(0 vi) component distribution 
(1032 A line) below z ^ 0.5, shown in the left panels of 
Figure [3] is the most commonly published low-z O VI ob- 
servable. It is relatively independent of spectral resolution 
and noise, so long as a given equivalent width systems is 
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Figure 3. Four of our models, d32n256vzwl50 with and without sub-SPH turbulence added, d32n256hzw075, and dl6n256vzwl50, are 
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Turbulence broadens the stronger Hi lines noticeably. 



identifiable. T08 performed a blind O VI search (i.e. search 
for the Ovi doublet alone without requiring corresponding 
Hi), explicitly differentiating between components and sys- 
tems. DS08 identify Ovi associated with Hi absorbers, and 
estimate they miss ~ 20% of O VI absorbers this way. Their 
larger sample has the advantage of reach ing down to 10 mA. 
This is especially important, because iTumlinson &: Fand 
<|2005l ) suggest that a turnover in this distribution indicates 
the extent to which metals are distributed to lower over- 
densities. The turnover in the DS08 data set is no t as pro- 
nounced as in the smaller iDanforth fc Shulll (|2005h dataset 
at 30 mA. We will argue that the declining numbers are 
weakly indicative of a metallicity gradient with density. Fig- 
ure[T]illustrates the decline in O VI at weak column densities; 
fewer regions exist where iV(Ovi) < 10 13 cm -2 (green re- 
gions in the bottom right) compared to Hi (bottom left), 
whic h is not observed t o turnover at iV(H i) < 10 13 cm -2 
(e.g. iPenton et ai1l2004l ). 

Our second observable is b(0 vi) as a function of 
N(Ovi). In the literature, this plot i s often used to d i agnos e 
the physical state of absorbers, as in iHeckman et~aH l|2002l) . 
That paper models O VI absorbers as radiatively cool- 
ing flows passing through the coronal temperature regime, 
finding a linear correlation between iV(Ovi) and 6(0 vi) 
corresponding to the cooling flow velocity, and indepen- 
dent of density and metallicity. They also show non-IGM 
Ovi absorbers (Ovi around local group gala xies and star- 
bursts ) follow this relation very closely, and iLehner et al.l 
(2006) finds that IGM absorbers agree with the same trend. 
However, this latter paper fails to match the predicted 
A r (Nevm)/A r (0 vi) ratios, and suggests th e IGM lines may 
instead be photo-ionized. IDanforth et alj (120061) a l so do 
not see O vi absorbers following the IHeckman et al.l (|2002l ) 
trend, but note that their lines may suffer from blending and 
exaggerated 6-parameters. 



We use only the highest resolution (STIS) data of in- 
tervening absorbers from T08, and plot these 55 absorbers 
spanning z = 0.13 — 0.50 in 0.2 dex bins in the central pan- 
els of Figure [4] with 1 a dispersions. We consider 6(0 vi) 
as a function of iV(0 Vi) (rather than vice-versa) because 
we will show that our absorbers are primarily photo-ionized 
with the underlying relation being density- 7V(0 vi). A clear 
upward trend is obvious in this data sample with lines 
iV(Ovi) > 10 14 cm~ 2 averaging 35.6 kms -1 . These lines 
cannot be explained by thermal broadening alone, because 
the thermal fc-parameter is only 17.7 kms -1 for the opti- 
mal O vi CIE temperature of » 300, 000 K. Hence there is 
extra broadening in these lines; possibilities include instru- 
mental broadening, Hubble flow broadening, and turbulent 
broadening, all of which we will explore. 

The final observable is N(0 vi) as a function of iV(Hl) 
for well-a ligned components , show n in the right panels of 
Figure H IDanforth fc Shulll l|2005l ) fit this relation with a 
iV(Hl) 0,1 power law, i.e. iV(Ovi) rises very mildly with 
iV(Hl). They thus argue that these two ions trace different 
IGM phases, with O vi tracing collisionally ionized gas, and 
consider this the multi-phase ratio. T08 also finds very little 
variation of iV(0 vi) with iV(H i), but show that single-phase 
photo-ionized models can reproduce much of this trend if 
metallicities are allowed to vary. They find their aligned ab- 
sorbers are well-described by the median of their robust 
Ovi sample, TV(Ovi) = 10 13,84 cm -2 , shown as a solid 
line. They further argue that well-aligned O vi with H I is 
most likely pho t o-ioniz ed; this conclusion is also reached by 
iThom fc Chenl (l2008al ). who argue that slightly increasing 
O VI strength with greater JV(H i) implies a declining ion- 
ization parameter with increasing gas column density. The 
specifics of the physics revealed by this trend are explored 
further in UXT\ 

We define "well-aligned" O vi and H I absorption as hav- 
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Figure 4. Nine models are fit to three Ovi observables: the cumulative EW distribution (left) with data from T08 (squares) and DS08 
(triangles), b(0 vi) as a function of N(Q vi) (c enter) over z = 0.15 — 0.5 from T08 components, and the JV(0 vi) as a function of A^Hl) 
(right) with a fit from lDanforth fc Shull j2005l) (dashed line) and the relation from T08 (where N(0 VI ) = 10 13 ' 84 cm -2 , solid line) along 
with two percentages indicating the fraction of well-aligned Ovi absorbers with JV(Ovi) = lo 13 - 5 " 140 (left) and JV(Ovi) = 10 140-150 
(right); the corresponding observed alignment fractions are 57 and 43% respectively. The top panels explore how different wind models fit. 
The middle panels explore other model variations such as a quasar-only ionization backg round scaled to th e H I ionization rate, a uniform 
Z = 0.1 Zq metallicity distribution, and non-equilibrium ionization ratios from lGnat fc Sternb erg (2007) for Ovi at CIE densities. The 
bottom panels show the addition of sub-SPH turbulence for all three wind models, which we suggest are small-scale velocities below our 
resolution limit in our simulations, yet provides the best fit to the data by broadening lines and increasing the EW of the strongest 
absorbers. Error bars show 1 cr errors in our sample on the left and right panels, while in the center panel they indicate 1 a dispersions 
to indicate the range of fe-parameters. Slight offsets in the x-axis values are for the ease of visibility. 



ing Sv = |f(Hl) — v(0 vi) | ^ 8 kms -1 , which corresponds to 
a Hubble-broadened spatial extent of 103 ± 3 kpc given the 
uncertainty in cosmological parameters from the WMAP5 
results. We compare our models to the T08 dataset by ap- 
plying the exact same criteria in searching for the closest 
Hi component to every Ovi component. We only consider 
absorbers at JV(Ovi) 10 13 ' 5 cm -2 , which is essentially a 



complete sample in S/N ^ 10 spectra. For the same reason, 
we require AT (Hi) ^ 10 12 ' 9 , so that an observed high-S/N 
sight line does not have a greater alignment fraction due to 
more detected weak Hi absorbers compared to a low-S/N 
sight line. A single H I component is allowed to be associated 
with multiple O vi absorbers; the converse is not applicable 
when we discuss the fraction of Hi absorbers aligned with 
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O VI later in l5.2l T08 applies an alignment criteria that con- 
siders instrumental resolution (see their §2.3.3), but we use 
a simple limit of 8 kms -1 because almost every component 
is determined to an error below this li mit. We use the T08 
dataset over the lThom fc Chenl \ 2008al ) dataset because the 
former is larger and spans a larger range of JV(Ovi). We 
make an exception for one strong system at z — 0.20266 in 
PKS 0312-770, since T08 does not fit Hi owing due to large 
uncertainties, and we i nstead use the parameters found by 
iThom fc Chenl (l2008bl ). This is an important system, be- 
cause it is one of the strongest Ovi absorbers with the 
strongest aligned Hi. 

In the upper left area of the right-hand panels in Fig- 
ure [4] we show two percentages, corresponding to the frac- 
tion of well-aligned Ovi components in two intermediate- 
strength 7V(0 vi) bins: 10 13 5 - 14 and lo 14 " 15 cm' 2 . For 
comparison, the T08 dataset has 57% and 43% of strong ab- 
sorbers well-aligned in the low and high bins respectively by 
our criteria. T08 considers 37% of their O VI systems aligned, 
however this is a very different consideration where such sys- 
tems have a single O VI component well-aligned with H I. We 
do not consider simple versus complex systems as defined by 
T08, because the definitions are not easily applicable to an 
automated procedure as we run on our simulations; compar- 
ing simulated and real systems is a challenge left for a future 
paper where the data reduction and analysis are replicated 
as closely as possible. 

In the following subsections, we step through each of 
our models, and assess how the variations in input physics 
impact observable properties. We emphasize the observables 
that differentiate amongst models in the discussion below. 
We save the interpretation of well-aligned absorbers for a 
detailed examination of ionization models in H4.3.1I 

3.2.1 Wind Velocity Strength 

The top panels of Figure [3] show a comparison between our 
three wind models, with our default specexbin parameters. 
For all three wind models, there is a clear dearth of high-IiW 
lines, and the observed line widths are progressively larger 
than predicted to higher N(0 vi). Hence wind strength does 
not affect the high-£W components nor their fc-parameters, 
which is an indication that the strongest lines are saturated 
and depend little on their actual metallicity as long as it 
is above a certain threshold. Since varying wind strength 
cannot explain these data, we require another explanation, 
which we will explore in the coming sections. ICen fc Fand" 
(2006) also find their simulations with and without super- 
winds make no difference at EW > 100 mA; however our 
SPH code versus their grid-based code treats the distribu- 
tion of metals on large scales very differently. In our case, 
Gadget-2 simulations without outflows barely produce any 
IGM absorbers (Springel & Hernquist 2003; OD06), thus we 
do not explore this case at low-z. 

Varying the strength of the winds has the largest ob- 
servable consequence for the number of low-iiW O vi com- 
ponents. Comparing vzw and Izw b est illustrates the sug- 
gestion of iTumlinson fc Fand l|2005h that the turnover seen 
at low EW is a consequence of declining IGM enrichment 
further from galaxies, and can b e used to const r ain th e ex- 
tent of the dispersal of metals. iFang fc Brvanl (|200ll ) find 
their simulations generate a turnover in the EW distribu- 
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Figure 5. The average metallicity per density for the three wind 
models in the top panel show that while all wind strengths form a 
metallicity-density gradient, stronger winds flatten this relation- 
ship. The volume filling factors at various limiting metallicities 
in the bottom panel indicate wind strength plays a major role. 
At Z = 0.1 Ziq, a factor of 27 reduction in volume filling factor 
across the three wind models is reflected in observations only by 
50% fewer weak components. 



tion at z — due to a metallicity-density gradient. Indeed, 
all our simulations produce metallicity-density gradient as 
we display in the top panel of Figure [5] OD06 also showed 
strong gradients at higher redshift. While the steeper gradi- 
ent of the Izw model yields an observational signature, the 
milder flattening of the hzw below overdensities of 10 has no 
observable effect. 

O VI does not distinguish the density-metallicity gradi- 
ent well; however, a key point is that our O vi absorbers arise 
from a fairly small volume filling factor of metals, and this is 
strongly governed by fj w i n d- By z = 0.25, our three simula- 
tions have volume filling factors where Z 0.01 Zq (deter- 
mined from oxygen) of 1.9%, 5.5%, and 18.0% as v W md in- 
creases (lzw~^vzw—>hzw). We determine these filling factors 
from lines of sight in physical space (i.e. no velocity smear- 
ing) using specexbin, plotting them in the bottom panel of 
Figure for various metallicities. Ovi generally arises from 
regions with Z ^ 0.1 Zq, where volume filling factors are 
only 0.3%, 1.3%, and 8.1% with increasing u w ind- The Izw 
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model differentiates between the amount of volume filled 
with half as many absorbers at EW = 1 5 mA and below, but 
the vzw and hzw models are indistinguishable despite even 
larger disparities in the volume filled. This disparity grows 
even further at the Zq limit; the hzw model enriches 18 x 
more volume than the vzw model, which fills only 1/15, 000th 
of a cosmological volume! These relatively small filling fac- 
tors in our simul ations are a consequen c e of t he fact that, 
as we showed in IDave fc Oppenheimerj (|2007T) and OD08, 
winds at late times enrich an increasingly smaller comoving 
volume, while many of the metals launched at early times 
reaccrete onto galaxies. Hence metals increasingly reside in 
compact high-overdensity regions at lower redshifts. 

It is remarkable that an order of magnitude change in 
the volume filling factor at the fiducial Ovi metallicity is 
observationally indistinguishable. It is not mysterious how- 
ever, because the mass-weighted oxygen abundance (only 
7% less in vzw compared to hzw) determines most trends 
in the observables. The Figure [1] slice of the d32n256vzwl50 
simulation illustrates that the volume enriched (upper right) 
does not correlate to where observable O VI is (bottom right; 
compare blue, pink, & white regions). 

The alignment of H I and O VI appears to provide no in- 
dication of the metallicity distribution either. The percent- 
age of aligned components is statistically indistinguishable 
in the two N(0 VI ) bins in the upper right panel of Fig- 
ure |3] as well, which could be considered surprising given 
how differently oxygen is distributed. We take this as an 
indication that the peculiar velocities are more important 
when considering aligned components rather than spatial 
distribution. It is disappointing that one cannot constrain 
the metallicity-density gradient or the volume filling factor 
from EW (Ovi), despite O VI being an ideal t racer of photo- 
ionized low-density metals (IDave et al.lll998h . 

3.2.2 Uniform Z = 0.1 Z Distribution 

In the middle row panels of Figure [4] we explore variations 
of post-run physics. Applying a uniform metal distribution 
(orange lines in the middle panels) goes in the opposite di- 
rection of how our models need to be modified to fit the EW 
distribution, since there are now more metals in low-density 
regions and fewer in high-density regions compared to our 
self-consistent enrichment models, which have metallicity- 
density gradients shown in Figure [S] Conversely, it illus- 
trates the possibility of explaining large 6-parameters with 
spatially extended structures undergoing Hubble expansion. 
Stronger components do show greater spatial broadening, 
with b ~ 30 kms" 1 at TV(Ovi) ^ 10 14 ' 1 cm -2 correspond- 
ing to spatial scales of ~ 400 kpc. With 100% volume filling 
factor of Z ^ 0.1 Z©, the number of 10 mA compo- 
nents only triples over the d32n256vzwl50 wind-distributed 
metals, which has a volume filling factor 80 x less; again, 
this shows that \ow-EW lines are a weak tracer of volume 
filling factor. Lines with EW ^ 100 mA are statistically in- 
distinguishable suggesting that these thin lines are already 
saturated with only 0.1 Zq. 

The slope of iV(Hi) — iV(0 vi) relation steepens due to 
more weak O VI absorbers, because there is an excess of such 
components near the O VI detection limit. This model does 
show a 17% higher incidence of well-aligned components 
with jV(Ovi) = lO 13 ' 5 " 14 cm -2 (78%) than any other 



model, which is a statistically meaningful difference. Met- 
als are smoothly distributed over the same fluctuations that 
create H I absorption leading to this higher fraction. The fact 
that the T08 dataset finds a lower fraction (57%) is an indi- 
cation that metals are distributed differently than baryons 
traced by Hi. In summary, applying a uniform metallicity 
does not appear to reconcile the models and observations. 

3.2.3 Ionization Background 

The middle row of Figure [3] illustrates the impact of varying 
the ionizing background. 37% fewer EW ^ 30 mA O VI lines 
create a more obvious turnover when the quasar-only back- 
ground is tripled (ibkgd-Q3), while leaving the Lya forest 
statistics nearly unchanged. The peak of photo-ionized O VI 
ionization fractions moves to higher overdensities due to the 
higher ionization parameter at 8.4 Rydbergs, resulting in 
O VI not tracing lower overdensities (~ 10) nearly as much, 
and leading to an underestimate of weak absorbers relative 
to the data. If future datasets from COS show a greater 
downturn in the EW distribution, then perhaps this harder 
ionization field is more appropriate. Little is know about the 
ionization field at 8.4 Ry dbergs, however the quasar+galaxy 
lHaardt fe~ Madau (2001) background works extremely well 
for both the Hi and Ovi observables we consider. 

The ibkgd-Ql background leaves the Ovi statistics un- 
altered, but creates significantly more H I absorption, in con- 
flict with observations. We do not plot this in Figure [4] as 
the only change is the H I — O VI relation moves rightward by 
nearly 0.5 dex due to greater Hi optical depths at a given 
overdensity. 

3.2.4 N on- Equilibrium Ionization 

Lastly, the middle row panels of Figure [4] illustrate t he im - 
pact of non-equilibrium cooling. iGnat fc Sternberg! (|2007f ) 
calculate non-equilibrium ionization states and cooling ef- 
ficiencies for rapidly cooling shock-heated gas without any 
ionization background. Since recombination lags cooling, the 
temperature range over which collisionally ionized O VI can 
exist extends below T — 10 s K. When we apply only their 
metallicity- dependent ionization fractions (their isochroic 
case) to where Ovi is primarily collisionally ionized (i.e. 
nh > 10~ 41 cm -3 at z — 0.25), we find a slight but im- 
measurable increase in O VI absorption from complex sys- 
tems tracing halo gas. Collisionally ionized O VI is primarily 
found near galactic halos, which have a small volume filling 
fact or at \ow-z (s ee i|5.3|) . 

ICen fc Fang] ||200d ) added non-equilibrium ionization 
directly to their simulations, and find an appreciable re- 
duction in Ovi tracing lower overdensities, because it is 
here that CIE timescales are longer than shock-heating 
timescales; Ovi remains in CIE longer in the equilib- 
rium case. Similar tw o-temperature behavior was shown by 
lYoshida et all (|2005l ) to have observational consequences in 
WHIM and ICM gas when applied to their cosmological sim- 
ulations. Hence it is possible that including non-equilibrium 
ionization self-consistently during the simulation may yield 
a larger difference. 

As an initial attempt towards this, we run the 
dl6nl28hzw075-noneq simulation with non-equilibrium 
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metal-line cooling from lGnat fc Sternberg! (|2007h during the 
run. However, this creates no appreciable difference in Ovi 
either. The broader metal-line cooling bumps actually ap- 
pear to make it slightly easier for baryons to cool into galax- 
ies, and there is a slight increase in Q*/Qb at z — (6.3% 
versus 5.7% from our dl6nl28hzw075 test simulation). How- 
ever, O VI observables are insensitive to such a small change 
in star formation, and the brown curves in the middle row 
panels of FigureUare indistinguishable from the assumption 
of CIE. 

Density-dependent non-equilibrium ionization effects 
with the inclusion of an ionization background could affect 
photo-ionized Ovi, especially given the fact that these met- 
als get into the diffuse IGM via shocked outflows. Recom- 
bination timescales exceed the Hubble time in such regions, 
and the two-temperature behavior of ions and electrons may 
have unforeseen consequences, possibly moving oxygen seen 
as O VI to higher ionization states, ft remains to be seen 
how big of an effect this is, and non-equilibrium effects on 
photo- ionized Ovi will be necessary to explore. 

Overall, non-equilibrium ionization effects may have 
some impact on weak Ovi systems, but is unlikely to pro- 
vide an appreciable change to {(-parameters as required to 
match the N(0 vi) — b(0 vi) relation. We plan to implement 
full non-equilibrium oxygen network during our simulation 
run to study this possibility further, but it seems unlikely 
to yield a large difference. 



3.2.5 Sub-SPH Turbulence 

We now consider the by-hand addition of sub-resolution tur- 
bulence, given by b tur b as described in Equations and [5] 
In this case, T08's observed N(0 vi) — 6(0 vi) relation in 
the bottom center panel of Figure [4] is reproduced by con- 
struction. However, an added consequence of the broader 
lines is that some of the lines that were saturated in the 
non-turbulent case now lie on the extended linear portion of 
the curve of growth. Their equivalent widths grow dramati- 
cally, and these models provide a good fit to the cumulative 
EW distribution at 50 mA and above. Our favored model, 
the d32n256vzwl50-bturb, matches the T08 EW distribu- 
tion as well as the 10 mA data of DS08. All the wind models 
shown in the bottom panels behave in the same way relative 
to each other as discussed in H3.2.1I 

A major caveat is that, as we show in £13.2.61 the over- 
all abundance of O vi absorbers is not completely converged 
with numerical resolution, in the sense that higher resolu- 
tion runs tend to show more absorbers at high EW. Hence 
the good agreement with d32n256vzwl50-bturb should be 
taken with some caution, and our wind model may need 
further refinement. The robust result here is that only the 
addition of turbulence is able to change the shape of the EW 
distribution to resemble that observed. 

We focus on the differences between turbulent and 
non-turbulent cases in Figure HJ] where we have included 
the d32n256vzwl50 and d32n256vzwl50-bturb models along 
with the higher resolution dl6n256vzwl50 model, which 
we discuss below. Turbulence blends components, primarily 
weaker ones in complex systems, which would be identified 
individually without turbulence; this is evidenced by 21% 



fewer 10-30 mA component^], and could be only a minor fac- 
tor in explaining the turnover at low-£W. When considering 
fi(Ovi) by summing all components between 10 12 ' 5 — 10 15 
using equation 6 of DS08, we find S7(0 vi) = 3.8 x 10~ 7 with 
turbulencfl only 5% more than without; turbulence should 
not affect column densities much, because the same amount 
of O vi remains in a sight line. 

Both simulations have the same amount of O VI along a 
line of sight, but Auto VP will often fit components with and 
without turbulence dramatically differently. An example is 
an absorber at z — 0.19 where the turbulently broadened 
line is fitted with iV(Ovi) = 10 14 ' 73 cm" 2 , b = 38 fans" 1 , 
and EW = 305 mA, versus iV(Ovi) = 10 14 ' 33 cm" 2 , 
b = 9.4 fans" 1 , and EW = 90 mA in the non-turbulent 
case. Without turbulence, much of the O VI remains hidden 
in the saturated line center. T08 finds no lines with at least 
10 14 ' 3 cm" 2 having 6-parameters under 33 kms" 1 using 
STIS. This is the typical absorber that greatly enhances the 
high end of the turbulent EW distribution, increasing it by 
2.8 x despite having the same underlying physical distribu- 
tion of O vi. along the sight line. 

Every single one of our mod els follows shows increa sing 
iV(Ovi) with AT (Hi), which as iThom fc Chenl (|2008al ) ar- 
gue indicates a decreasing ionization parameter for stronger 
Ovi, and suggests a photo-ionized origin. The aligned ab- 
sorber percentages in the lower bin are only 2% apart 
with and without turbulence, and agree statistically with 
the T08 percentages. There are more aligned absorbers at 
iV(Hl) > 10 14 ' 5 cm" 2 , because of the blending of separate 
Hi components in very dense regions due to added turbu- 
lence. 

It is worthwhile to ask whether it is fair to add such a 
"fudge factor" in our models just to fit the data. In its favor, 
it does separately alter the shape of the EW distribution as 
desired, and there seem to be no other effects capable of 
doing so. Still, it is an unsatisfyingly ad hoc addition, and 
in t|4.4l we attempt to further explore the physics responsible 
for bturb, and understand their implications for small-scale 
motions in the IGM. 



3.2.6 Resolution Convergence 

We use the dl6n256vzwl50 simulation only to test resolu- 
tion convergence, and not to explore the general proper- 
ties of Ovi absorbers, since it does not contain a statis- 
tical sample of galaxies at M* nor a representative vol- 
ume of low-z environments. Although the dl6n256vzwl50 
simulation produces 40% more star formation, this is not 
evidenced in the low end of the EW distribution in Fig- 
ure [6] where there exist 26% fewer absorbers between 10- 
50 mA than the d32n256vzwl50 simulation. The excessive 
late-time star formation in the dl6n256vzwl50 simulation 
(as mentioned in §2.1\ enriches mainly higher overdensities 
probed by stronger Ovi absorbers (50% more 100 mA Ovi 



1 The cumulative EW distribution only shows a 7% difference 
summing all absorbers ^ 10 mA. 

2 This value compares favorably with DS08's value of 4.1 ± 0.5 X 
10" 7 above a similar limit (30 mA), but we do not wish to draw 
too much into this comparison, since our simulation resolution 
limits may be hiding more O VI, see £|3.2.6I 
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Figure 6. Three models with vzw winds are plotted against the same three Ovi observables in Figure [4] The d32n256vzwl50 and 
d32n256vzwl50-bturb are plotted on the same panels to directly show how turbulence affects observational signatures, while the 
dl6n256vzwl50 model explores the effect of increasing mass resolution by 8x. 



absorbers) and lower ionization metal species. Although the 
dl6n256vzwl50 simulation is resolving more galaxies down 
to ~ 10 s Mq, these galaxies do not seem to be a signif- 
icant contributor to the IGM Ovi. Instead, the difference 
for \ow-EW absorbers appears to depend on the hydrody- 
namics of outflows at different resolutions- winds from the 
same mass galaxy at higher resolution do not enrich as effec- 
tively low density regions where weak Ovi absorbers arise. 
This could be a result of more resolved substructure lead- 
ing to greater hydrodynamic slowing at higher resolutions. 
We caution that the differences could also be partly numeri- 
cal artifacts inherent to SPH, such as excessive viscosity and 
poor two-phase medium separation. Furthermore, there is no 
evidence we can resolve the sub-SPH turbulence any better 
at x8 greater resolution as evidenced by the fa-parameters 
showing nearly no change between the two simulations; this 
turbulence must be added explicitly at cosmological resolu- 
tions. 

The N(Ul)/N(Ovi) ratio in right panel of Figured] 
again shows essentially no difference between the two reso- 
lutions. There are more strong aligned Ovi absorbers above 
iV(Ovi) = 10 14 ' 1 cm -2 due to greater enrichment of over- 
dense regions in the dl6n256vzwl50 box. 

We do not explore our d64n256vzwl50 box, because (1) 
this simulation significantly underestimates SF at z > 2 as it 
cannot resolve below M* at these redshift, and (2) the mean 
interparticle spacing is greater than the sizes of absorbers 
(50 — 100 kpc) we find in the next section. The 32/i _1 Mpc 
256 3 is above the critical resolution where we can resolve 
the production of metals in O VI absorbers and the sizes we 
derive for these absorbers. 

Overall the resolution convergence of O VI absorber 
properties, while mostly statistically consistent in Figure 
[6] indicate some numerical issues within our simulations. 
It appears unlikely that modest increases in resolution will 
dramatically change the basic conclusion that sub-SPH tur- 
bulence is needed to match O VI observations. Running sim- 
ulations with sufficient resolution to study small-scale turbu- 
lence while properly modeling nonlinear structures to z = 
is well beyond present computing abilities. We further stress 
the need to run galactic-scale simulations in the future to 



fully resolve the absorbing structures within haloes and in- 
tragroup gas where the strongest O VI absorbers appear arise 
in our simulations. 

Finally, concerning box size convergence, it is true that 
the 32ft _1 Mpcs does not contain larger modes that could be- 
come non-linear by z — 0; however we will show that we are 
primarily exploring the diffuse IGM around sub-M* galax- 
ies, which are statistically well-sampled in our simulation. 
This is evidenced by the fact that our O VI and H I statistics 
in the dl6nl28vzwl50 box are statistically consistent with 
those of the d32n256vzw!50 box. 



4 PHYSICAL CONDITIONS OF Ovi 
ABSORBERS 

From here on we focus on our favored model with sub- 
resolution turbulence added, d32n256vzwl50-bturb, and ex- 
amine the physical nature of O vi absorbers. We first exam- 
ine the density, temperature, metallicity, and absorber size of 
IGM Ovi components, finding that most observed Ovi sys- 
tems are consistent with being photo-ionized. Shorter cool- 
ing times of a clumpy metal distribution are argued as the 
driver for Ovi to reach such temperatures. We then dis- 
cuss the correspondence of H I with O VI, using the N(R i) — 
iV(Ovi) relation to explore various photo-ionization, colli- 
sional ionization, and multi-phase ionization models. Lastly, 
we justify adding turbulent velocities on physical grounds, 
explaining why it should not be surprising that turbulence 
provides a dominant component to many IGM absorber pro- 
files. This will lead us to the next section where we examine 
the location and origin of O VI relative to galaxies. 

Physical parameters ascribed to absorbers such as den- 
sity and temperature can be weighted in multiple ways. Pos- 
sibilities include H I and O VI absorption weighting, or just 
normal mass weighting. We weight such quantities using the 
formula 



X{W) 



Y^wxx 



(7) 



where X is the property and W is the weight. For exam- 
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pie, if a quantity is weighted by Ovi absorption, W = 
m x Zo x /(Ovi), where m is mass of the element (e.g. 
SPH particle), Zo is the mass fraction in oxygen, and 
/(Ovi) = Ovi/O (the Ovi ionization fraction). X is of- 
ten temperature or density itself. The difference between 
Hi and Ovi- weighted quantities is critical for understand- 
ing the multi-phase nature of absorbers. 

We examine 70 continuous lines of sight with S/N — 20 
between z = 0.5 — > 0; we call this is our high qual- 
ity sample. 1511 Ovi components are identified above 
iV(Ovi) = 10 12 - 5 cm -2 . This sample allows us to obtain 
reasonable statistics of all absorbers between iV(Ovi) = 
10 12 ' 9-14 ' 9 cm" 2 , which we will refer to as the "observed 
range" since nearly all observed IGM O VI absorbers lie in 
this interval. 



4.1 A photo-ionized model for Ovi absorbers 

We examine the density, temperature, metallicity, and sizes 
of Ovi absorbers based on the observational quantity 
TV(Ovi). Figure [7] (top two panels) plots the median Ovi- 
weighted density and temperatures at the line centers. Over 
two decades of column density covering the observed range, 
Ovi absorbers show a steady increase in tih, while the tem- 
perature stays nearly constant at T ~ 10 4 ' 15 K. Also plotted 
in dashed lines are the 16th and 84th percentile values cor- 
responding to the 1 a dispersions. The scatter in density 
is significantly larger, a = 0.3 — 0.4 dex, compared to the 
scatter in temperature, < 0.2 dex over much of the observed 
range. Our first main conclusion is that the vast majority of 
O VI absorbers at all but the very strongest column densities 
are photo-ionized, with temperatures less than 30,000 K. 

We already presented the fit to the density as a function 
of JV(Ovi) in Equation 31 which fits the relation quite well 
regardless of added turbulence (magenta dashed line). The 
increase in density alone can explain 66% of the O VI column 
density increase, with all else being equal. The remaining 
increase must arise from variations in metallicity, ionization 
fraction, and absorber length. Metallicity is expected to ex- 
plain some of the increase, because our simulations show a 
consistent metallicity-density gradient, as shown in Figure[S] 
Fitting the points in the third panel of Figure [7] yields the 
relation 



log[Z ] = 0.45 x log[iV(Ovi)] 
for AT(Ovi) = lo 12 - 9 " 14 - 9 



,96, 



(8) 



cm - . Although density (66%) 
and metallicity (45%) variations can account for all of the 
dependence of N(0 vi), we will show that ionization fraction 
and absorber length also vary. 

The metallicities reported here always exceed the fidu- 
cial value of 0.1 Zg often assumed when modeling low- z O VI 
absorbers. Our absorbers are generally found to lie between 
0.15-1.0 Zq, which is consistent wi th metallicity determi- 
nations from aligned H I absorbers (IProchaska et al.l 1 20041; 
iLehner et al. 2006; ICooksev et al. I l2008l. ; T08). We contrast 
the oxygen- weighted metallicity of Ovi absorbers (Equa- 
tion [8]) to the mass- weighted oxygen metallicity-density gra- 
dient averaged over the entire simulation box at z — 0.25: 



Iog[Zo] = -4.26 + 0.361og[p/p] + 0.075(log[p/p]) 2 



(9) 
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Figure 7. Physical properties as functions of O VI column density 
in the d32n256vzwl50-bturb model are plotted a solid black lines 
with points in the top 3 panels for all our simulated absorbers in 
70 high-quality lines of sight below z = 0.5. Dashed lines corre- 
spond to lcr dispersions. The long-dashed magenta line is the fit 
from Equation [4] We also plot the simulations-average metallic- 
ity as a function density (converted to TV(Ovi) via Equation [4] 
purple long-dashed line) and the metallicity one would infer us- 
ing aligned iV(H 1) — iV(0 vi) absorbers as in Equation [15] (orange 
dashed-dot line). We consider Ovi fractions from CLOUDY at 
three photo-ionized temperatures (2nd to last panel), which sug- 
gest photo-ionized Ovi absorbers are typically 50-100 kpc long 
(last panel). Ovi in CIE will produce more compact absorbers at 
higher JV(Ovi) (red lines), which appear to trace gas near and 
within halos. The Hi-derived absorber length (orange dashed-dot 
line) assuming O vi traces the same gas as H I appears to be a 
poor approximation for strong absorbers, suggesting such aligned 
absorbers are multi-density. 



for p/p = 1 — 3000, which is a fit to the green curve in the 
top panel of Figure [5] This shows a 34% rise over the same 
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N(0 Vl) range, and is represented by the purple dashed line 
in the middle panel of Figure [7] The metallicities of O VI 
absorbers are on average 4 — 6x higher than the general 
IGM over the observed range. This is a clear sign that Ovi 
absorbers trace inhomogeneously distributed metals. The 
low scatter in the oxygen metallicity-density relation indi- 
cates that O VI absorbers rarely trace gas with average IGM 
metallicity. Hence the global metallicity-density relationship 
is a poor descriptor of the enrichment level in metal-line ab- 
sorption systems. 

The fourth panel in Figure [7] shows the O VI ioniza- 
tion fraction, computed by assuming the density— iV(0 vi) 
relation shown in th e top two panels , for various tem- 
peratures, assuming a lHaardt fc Madaul (|200lf ) quasar-only 
background. For a typical Ovi temperature of T — 10 4 ' 2 
K, /(Ovi) changes significantly from 0.23 to 0.029 over the 
observed range. It depends little on temperature for lines 
above 10 14 cm -2 . The peak ionization fraction for colli- 
sional ionization (T = 10 5 ' 45 ) is around 0.2, but it is only 
approached for the strongest lines. This is another reason 
why over most of the observed range, photo-ionization is 
the dominant mechanism for O VI absorbers. 

Thus far, the dependence of iV(0 vi) using photo- 
ionization is partially explained by density (66%) and metal- 
licity (45%), but the declining ionization fractions above 



model. In other words, the photo-ionized nature of IGM O vi 
absorbers is a highly robust prediction of our simulations. 



4.2 Cooling Times 

The inhomogeneous nature of the metal distribution is vi- 
tal for the photo-ionizational explanation of Ovi, because 
metal-line cooling (which we implemented in Gadget-2 in 
OD06) is much more effective when metals are concentrated. 
Consider the cooling time of gas in CIE at z = 0.25: 



r cool = 4.12 x 10" 



T 



5A(T, Z) 



yrs 



(11) 



10" 



cm anti-correlate with N(0 vi) (-45%). The larlyj a strong system w ith N(0 vi) = 10 1 



remaining dependence (34%) therefore must be in the ab- 
sorber size (£(Ovi)). Absorber size is calculated in the bot- 
tom panel of Figure[7Jby dividing the column density by the 
number density of O vi: 



1(0 vi) = 



JV(Ovi) 



n H (0 vi) 



,/h 



2a /(Ovi) 

mo K ' 



(10) 



where the /h = 0.76 is the mass fraction in hydrogen, and 
is the ratio of the atomic weights (0.0625). Absorbers 
remain between 50-100 kpc at AT(Ovi) < 10 14 cm -2 where 
/(Ovi) w 0.20 — 0.25, and rise to many hundreds of kpc 



at AT(Ovi) 



10 1 



cm assuming photo-ionization. The 



large required path length for the very strongest absorbers 
makes photo-ionized absorbers rarer, and makes collisionally 
ionized systems (T = 10 5 ' 45 ) with smaller sizes (red curve) 
more viable. We will discuss such systems further in £15.31 

In summary, virtually all O vi IGM absorbers are photo- 
ionized. Weak absorbers trace filaments with overdensities 
in the range of 10 — 20; intermediate absorbers trace over- 
densities 20 — 100; and strong absorbers trace gas within or 
near galactic halos. We stress that some collisionally ion- 
ized O vi absorbers do exist, particularly for the strongest 
absorbers where a much shorter pathlength can yield suffi- 
cient iV(Ovi), but these are better viewed as galactic halo 
rather than IGM absorbers. Strong O vi absorbers associ- 
ated with the MW halo indicate collisionally ionized oxygen, 
because otherwise their long pathlengths w ould be unable 
to fit within the halo l|Sembach et al. |2003l ); our absorber 
length analysis agrees with this assessment. Since our sight 
lines occasionally pass through such halos, we discuss such 
absorbers in £|5 . 3 1 when we consider environment. Finally, 
while we have used the d32n256vzwl50-bturb model, the 
trends identified here are nearly identical without sub-SPH 
turbulence added, and are essentially independent of wind 



where A is the cooling rate (ergs s _1 cm 3 Let us consider 
the cooling time of a typical 10 13 cm -2 absorber (corre- 
sponding to overdensity, S — 13), assuming it cools from 
3 x 10 5 K. If we use Z from Equation [5] (i.e. we assume the 
mean IGM metallicity) and assume solar abundance ratios 
for all species, the cooling time is 42 Gyr. Conversely, using 
the metallicity of O vi absorbers from Equation [5] yields a 
cooling time of 10 Gyr. Hence weak absorbers cannot cool 
significantly within a Hubble time at the mean metallicity, 
but can do so at the mean Ovi absorber metallicity. Simi- 

m" 2 (<5 ~ 100) 

has cooling times for a mean IGM metallicity and for Ovi 
absorber metallicity of t coo i = 2.9 and 0.6 Gyr, respectively. 
Of course, densities in the early Universe were higher by 
(1 + z) 3 , and cooling times depend inversely, but metallici- 
ties were lower, so these cooling times are probably overes- 
timates. Still, for observed O vi absorbers the cooling time 
is generally shorter than a Hubble time (rnubbie) when the 
Ovi absorber metallicity is considered. 

Weaker absorbers take longer to cool, but as we will 
show in H5. II weak absorbers trace some of the oldest met- 
als injected at high-z, while stronger absorbers trace more 
recent metal injection. This means that all but the very 
strongest O vi absorbers have had plenty of time to cool 
to photo-ionized temperatures since their injection into the 
IGM. Physically, this is why the clumpy metal distribution 
where Z ~ 5Z, causes most Ovi absorbers to be photo- 
ionized. 

In practice, superwind feedback shock-heats metals to 
higher temperature (T > 10 s K) in denser environs than 
those traced by photo-ionized Ovi. The result is a range 
of cooling times with lower metallicities more likely to have 
Tcooi > THubtie and remain in the WHIM or hot IGM, and 
higher metallicities often with Tcooi THubble- This lat- 
ter case includes metals traced by photo-ionized O VI. The 
strength of the feedback is a sensitive determinant for which 
(p, T) path an SPH particle will follow, with stronger feed- 
back putting many more metals into the WHIM and hot 
IGM (OD08). Our main point here is to show that metals 
heated to WHIM temperatures at O vi absorber metallicities 
will usually cool to photo-ionized temperatures by today. 



3 See Gnat & Sternberg (2007) for a derivation of the general 
form of this formula. 
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Figure 8. The alignment fraction of Ovi absorbers with Hi at 
various separations in two bins covering intermediate absorbers 
(JV(OVI) = lO 13 - 5 " 14 and io 14 - " 14 - 5 cm" 2 ). The T08 data is 
compared to the S/N = 10 sight lines of the d32n256vzwl50- 
bturb simulation, showing nearly perfect agreement in the lower 
bin (top) and adequate agreement in the upper bin (bottom). The 
high aligned absorber fractions are an indication of photo-ionized 
Ovi. 

4.3 Alignment with Hi 

We have thus far discussed O VI absorbers independently 
of Hi. Any model for Ovi absorbers should also reproduce 
the general properties of aligned Hi absorbers, in addition 
to reproducing global Lycv forest properties (e.g. IDave et al.l 
1 19991 ; IDave fe Trippll200ll ). We start by emphasizing that 
a vast majority of O VI absorbers are aligned at some level 
with Hi. However, from the Hi- and Ovi-weighted density 
and temperature of well-aligned absorbers, we will demon- 
strate that aligned absorbers often arise in different phases 
of gas. We introduce the idea of multi-phase photo-ionized 
models with different density and temperatures for the O VI 
and Hi components. We then estimate Hi-derived absorber 
length, showing these can differ from the O VI lengths often 
because they are tracing different phases of gas. Lastly, we 
fit multi-phase ionization models to the iV(Hl) — iV(Ovi) 
plane, to understand what types of absorbers populate vari- 
ous regions. This can aid in interpretation of data sets such 
as T08 and upcoming data from COS. 

The alignment fraction as a function of velocity sep- 
aration for this model, as for most of our models, agrees 
well with observations, as we show in Figure [8] for interme- 
diate absorbers split into two bins. Well-aligned absorbers 
(8v ^ 8 kms -1 or 100 kpc in a Hubble flow) comprise 57% 
of all Ovi absorbers in the d32n256vzwl50-bturb simula- 



tion when considering intermediate absorbers, in statistical 
agreement with T08 (54%) (see §3/2$. We use the S/N = 10 
lines of sight in this case, because it is more similar to the 
data S/N; the high quality sample has a 62% alignment frac- 
tion, since more components are identified. The agreement 
in the iV(Ovi) = io 13 - 5 " 14 cm -2 bin is nearly perfect. 
There are fewer aligned absorbers in the data for the higher 
column density bin (TV(Ovi) = io 14 " 14 - 5 cm" 2 ), but this 
is not a statistically significant difference since the observed 
sample has only 14 absorbers. 

Misaligned absorbers are an important population, be- 
ing a signature of collisionally ionized O VI as we discuss 
in £15.31 Our discussion in that section suggests strong ab- 
sorbers have a greater chance of being collisionally ionized; 
however we argue in the following subsections that the high 
number of aligned absorbers in both the data and simula- 
tions is a strong indication of photo-ionized O VI. 

4-3.1 Density and Temperature 

While we have argued that the vas t majority of Ovi ab - 
sorbers are aligned at some level (e.g. lThom fe Chenll2008al ) ■ 
do such aligned absorbers trace the same underlying gas? To 
answer this we plot the H I and O VI- weighted densities and 
temperatures for well-aligned absorbers against each other 
in Figure [9] H I and O VI trace the same density gas be- 
low 7ih ~ 10 -4 ' 8 cm -3 , but above this the absorbers tend 
to become multi-phase in density. The declining ionization 
fraction of photo-ionized O VI makes it a poor tracer of gas 
inside halos at overdensities above a couple hundred. These 
rare instances are less represented in well-aligned Ovi ab- 
sorbers due to the large peculiar velocities of various compo- 
nents. However, the spatial extent of regions at overdensities 
~ 500 traced by strong Hi is small (< 100 kpc), and gas at 
lower overdensities tracing strong Ovi absorption is often 
close by and therefore appears aligned. 

The weighted temperatures show significantly less align- 
ment, and reveal something that may be counter-intuitive 
- metal-enriched gas is cooler than unenriched gas. It is 
counter-intuitive in the sense that to reach the IGM via su- 
perwind feedback, gas is driven at high velocities and there- 
fore must shock heat. However, stronger metal-line cooling 
outweighs this by low-z. The Hi temperatures span a wider 
range, just falling short of the T = 10 K limit, which enters 
the realm of BLAs (fe(Hl) = 40 kms -1 ). The two tempera- 
tures show that even in aligned absorbers at similar densi- 
ties, Ovi and Hi are not tracing the same gas. Again, we 
suggest that O VI is arising from a clumpy distribution that 
picks out enriched regions of the more smoothly varying Lycv 
forest. Several rare instances of O VI in CIE is aligned with 
cooler gas with Hi absorption. 

4.3.2 Absorber Size 

Before introducing ionization models predicting line ratios, 
we must consider the Hi-weighted absorber size, which is 
different than /(Ovi). We show the Hi- weighted density for 
all Hi absorbers, even those without Ovi, in the top panel 
of Figure [10] A fit to the median of these absorbers between 
N(Ul) = io 12 - 14 5 cm" 2 gives the relation 

log[n H ( cm -3 )] = -15.68 + 0.77 log[JV(Hi)] (12) 
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Figure 9. The multi-phase nature of well-aligned absorbers is shown by plotting Hi and Ovi-weighted temperatures and densities 
against each other (left & right panels respectively) from 70 sight lines. The column density of absorbers are indicated by both the size 
and color of the error bar with 7V(Hl) in the x direction and ./V(Ovi) in the y direction. Densities diverge above = W~ 4 8 cm -3 
where O VI more often traces gas outside halos while H I traces gas within. The dashed line indicates the relation we use between the two 
densities. Temperatures are generally cooler for metal-enriched gas indicating metal-line cooling plays an important factor in gas, which 
indicates an inhomogeneous distribution of metals. 



at (z) = 0.25, and well describes the lower density bound- 
ary of Hi absorbers over 7 or ders of mag n itude . This re- 
lationship is also derived by iDave et al.l (|l999h . finding 
log(n H ) = -15.1 3 + 0.7 log[7V(Hi)] at z = 0.25, and analyt- 
ically by ISchavd (|200ll ) by assuming the characteristic size 
of an absorber is well-described by the IGM Jeans length, 
log(n H ) = -14.68 + 2/3 logfiVfHl)]. Our current simula - 
tions (which are much improved over IDave | Tripdl200ir i 
find an overall slightly steeper dependence and a higher nor- 
malized «H. The upward scatter in the top panel of Figure 
1101 are the absorbers likely to be part of com p lex H I sys- 
tems; similar scatter is seen by IDave fc Tripd l|200ll ). The 
well-aligned absorbers plotted in the bottom of Figure \W\ 
are consistent with equation[T2](red line), showing that they 
generally are not a biased population of Hi absorbers, ex- 
cept that stronger H I is more likely to have aligned O vQ- 

We take Equation 1121 as the relation between iV(Hl) — 
«h for all aligned absorbers, despite the scatter, and use 
the Hi fraction for gas at T = 10 4 ' 3 K photo-ionized by 
a quasar-dominated ionization background, log[/(Hl)] = 
0.17 + 0.99 log (nii), to get the neutral hydrogen column 
density relationship: 



log[7V(H) ( cm" 2 )] = 20.19 + 0.31 log(n H ). 



(13) 



Dividing iV(H) by nu yields the absorber length, 

log[Z(Hl) (cm)] = 20.19 - 0.69 log(n H ). (14) 

We plot i(Hl), using the assumption that nH(Ovi) = 
wh(Hi) as an orange dash-dot line in Figure [7] to con- 
trast with the O vi-derived absorber size, 50 — 100 kpc at 
AT(Ovi) < 10 14 ' 5 cm" 2 . Below JV(Ovi) = 10 140 cm" 2 , 
it seems fair to assume the same densities since nn < 



4 See the end of £|5.2I for the predicted fraction of H I absorbers 
with Ovi. 



10 -4 ' 5 cm -3 , and the resulting Ovi absorber lengths fit 
within the Hi absorber lengths. Above 10 140 cm -2 , i(Hl) 
becomes shorter than i(Ovi) for several reasons. First and 
perhaps most important, the densities where the two species 
arise diverge above nu = 10~ 4 ' 8 cm -3 , with Ovi arising 
from lower density gas distributed over larger lengths. Sec- 
ond, stronger absorbers are less likely to be aligned due to 
larger peculiar velocities, making this comparison less mean- 
ingful. Third, strong absorbers are more likely to have met- 
als at CIE temperatures, raising the ionization fraction and 
lowering i(Ovi). Finally, the successive scatter in each step 
above to get i(Hl) adds up when considering the spread in 
temperatures. 

As an interesting aside, we investigate the "iV(H) 
conspiracy" for low-z metal line systems presented by 
IProchaska et alj <|2004T ) ■ They find six O VI absorbers in the 
PKS 0405-123 spectrum averaging N(H) = 10 18 ' 7 cm -2 
with a scatter of only 0.3 dex. In fact, this is straight- 
forwardly understood in the context of our photo-ionized 
model for O VI absorbers. The iV(H) conspiracy comes about 
because O VI arises mostly from photo-ionized overdensi- 
ties (~ 10 — 200 at z = 0.25), which leads to a similarly 



small range, N(H) 



10 



18.50-18.91 



cm (using eq. [IB")) , 



in our model. This r ange is in excellent agreement with 
IProchaska et al.l (|2004T ) , providing further evidence that O VI 
in our simulations is tracing the cosmic web as observed. 
iV(H) column densities below 10 18 ' 5 c m" 2 generally indica te 
photo- ionized absorbers according to lRichter et al.l (|2006T ) . 

A common method to estimate oxygen metallicities of 
aligned absorbers from observations is 



JV(Ovi) /(Hi) mo 

Zo ~ IvThiT x 7(o^) x ^ x /H - 



(15) 



Unfortunately, this may significantly mis-estimate the true 
metallicity. To illustrate this, we plot this "observed metal- 
licity" as the orange dash-dotted line in the central panel 
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Figure 10. The Hi column density is related to njj over 7 
decades as shown by the fit to the median density between 
AT (Hi) = 10 120 ~ 14 ' 5 cm" 2 . Hlabsorbers aligned with O VI follow 
this relation as well, however in both cases scatter exists as weak 
absorbers sometimes trace high densities in complex systems. The 
Tin — N(R i) relation is used to determine the H 1-derived absorber 
size. 

of Figure [7] At weak Ovi column densities, the metallicity 
is underestimated, while at above iV(Hl) > 10 14 cm -2 it 
overestimates the actual metallicity. 

In summary, the H I properties of well-aligned O VI ab- 
sorbers are fairly similar to general Hi absorber properties. 
H I and O VI absorber sizes can be substantially different in 
stronger absorbers, because even aligned H I and O VI arise 
in different phases of gas. Finally, A/(H) remaining relatively 
invariant over a large range of A/ (Ovi) is expected if well- 
aligned absorbers arise from typical H I absorbers. 

4-3.3 Ionization Models 

Having the length scale of aligned absorbers, along with den- 
sities and temperatures, we now have everything necessary 
to fit ionization models to O VI — H I line ratios. This is por- 
trayed in Figure [TT1 using the JV(Hi) — iV(Ovi) multi-phase 



plane considered bv lDanforth fc Shulll (|2005h and T08. Note 
that these authors plot JV(Hi) versus A/(Hi)/AT(Ovi); we 
prefer to plot the column densities against each other di- 
rectly and not the ratio, because it is simpler and makes 
trends clearer. Small points represent simulated aligned ab- 
sorbers. We plot colored lines for various two-temperature, 
multi-phase models where oxygen is either photo-ionized 
(T(Ovi) = 10 4 - 2 K) or in CIE (T(Ovi) = 10 5 45 K), and 
Hi is one of several different temperatures. Our models 
use the z — 0.25 metallicity-density relation (Equation [9]l 
shown on top (dashed lines), but we also consider metallic- 
ities 5xZ (solid lines) as more relevant given the Ovi ab- 
sorber metallicities from ij4.ll It is difficult to create photo- 
ionized Ovi absorbers with column densities much greater 
than 10 14 ' 3 cm~ 2 using the 5x2 relation, because /(Ovi) 
declines as AT (Hi) rises creating a maximum in Af(Ovi). 

The CIE models (or ange, red, fc mage n ta) fo llow the 
behavior highlighted by iDanforth fc Shulll (|2005l ) where 
N(Ovi) is independent of A^Hl), which sweeps across a 
broad swath of this plot. In our case, the Hi tempera- 
ture ranging from T = 10 4 ' 2 — > 10 5 ' 5 K is the variation 
that allows this scenario to explain most of the T08 well- 
aligned absorbers (Sv ^ 8 kms -1 separation, blue squares) 
while O VI is always collisionally ionized. This scenario seems 
somewhat implausible at face value. Meanwhile, the photo- 
ionized models with two choices of T(Hl) (cyan and green) 
bracket the majority of H i-derived temperatures in Figure[7|] 
and also can explain most, but not all, of the T08 absorbers. 
Lastly, we introduce a multi-density model assuming T = 
10 4 ' 2 K (purple) for densities above nn = 10~ 4,5 cm -3 
where we assume different absorber lengths for O VI and 
Hi by plugging in the divergent density relationship from 
Figure (dashed line) into Equation 1141 This results in 
longer absorber pathlengths for Ovi than Hi. Unless oth- 
erwise noted, we assume aligned absorbers are described by 
photo- ionization. 

We compare the number densities of O VI systems to 
aligned T08 absorbers in three bins denoted by the vertical 
gray lines in Figure 0- AT(Hi) < 10 13 4 , 10 13 ' 4 " 15 , and 
> io 15 ' 0-17 ' cm -2 , and we overlay our high quality sample 
of 70 lines of sight (black dots, T(Ovi) < 10 5 K; red dots, 
T(Ovi) > 10 5 K). Our sample covers Az = 35 with ~ 
90% completeness at TV(Ovi) = 10 13 cm -2 . We directly 
compare with the T08 dataset only for absorbers of at least 
intermediate strength, for which they have Az = 2.62 (their 
effective pathlength at EW = 30mA). 

Beginning with the strongest Hi bin first, T08 finds 
5 such absorbers (1.9 Az -1 ), which is nearly 2 a more 
than we find (14, 0.4 Az -1 ). These absorbers correspond to 
gas in and around halos, and rarely fall below N(0 vi) = 
10 13 ' 6 cm -2 in either set. Two of T08's absorbers and 
many of our absorbers are above the solid T(Hl) = 10 4 ' 2 
K photo-ionized model (cyan) which seemingly rules out 
photo-ionization, because higher metallicities are required. 
We suggest that a subset of these absorbers are best de- 
scribed by the multi-density photo-ionized model (purple), 
with O VI arising from longer absorbers outside a halo, while 
H I traces denser halo gas. The higher observed frequency of 
strong absorbers in the data may indicate non-uniformity 
in the ionization background, where halo substructure unre- 
solved in our simulations increases the number of absorbers 
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Figure 11. The ./V(H i)— 7V(0 vi) plane shows the aligned absorbers from our high-quality lines of sight (black points, T(0 vi) < 10 K; red 
points, T(0 vi) ^ 10 5 K) along with absorbers observed by T08 (open blue squares). Several ionization models (colored lines) are shown 
using either the metallicity-density relationship. Z, in our z = 0.25 d32n256vzwl50 simulation (dashed lines) or 5xZ (solid lines) since 
metals in Ovi absorbers appear to be significantly above the cosmic average. The top panel shows this relationship translated to ./V(Hl) 
using Equation [12] and the relation n H = 3.63 X 10~ 7 p/p cm" 3 , which applies at z = 0.25. Collisionally ionized O VI (T(0 vi) = 10 5 45 K) 
aligned with Hi at various temperatures (orange, red, &; magenta) can explain most absorbers, but the absorbers are vastly photo-ionized 
and better explained by cyan and green ionization models (T(0 vi) = 10 4 ' 2 K). The purple ionization models are multi-density with Ovi 
arising from lower densities than Hi when njj 10 — cm -3 . 



with iV(Hl) > 10 15 cm 2 , which our simulations under- 
predict. 

It is significant that few of our simulated absorbers and 
none of T08's correspond to values lie below Z in either of 
the single-phase models (green and cyan dashed lines) above 



iV(Hi) = 10 15 cm 2 , although these shoul d be detectable 
in the data; this behavior also exists in the iThom fc Chenl 
(2008a) dataset. The Hi often traces halo gas while the 
photo- ionized Ovi arises from a longer pathlength of gas 
at 8 ~ 100, increasing the chance of intersecting with a 



20 B. D. Oppenheimer & R. Dave 



clump of metals well above Z. Only 29% of our strong H I 
absorbers have well-aligned O VI, but 87% of them have O VI 
within 80 kms -1 (see the end of £15.21 for more statistics). 
Our interpretation is there is almost always O VI with strong 
Hi, but the Ovi is much more extended and peculiar veloc- 
ities dominate, resulting in well-aligned strong Hi absorbers 
often being chance alignments. 

We find a vast majority of aligned absorbers in the in- 
termediate bin (8.5 Az _1 ), within the Poisson errors of T08 
(7.3 Az' 1 ). Most of our absorbers including those down to 
our detection limit fall between 1 — 10Z, which should hap- 
pen according to our photo-ionized O VI model. These ab- 
sorbers tracing 5 = 10 — 100 are more likely to be associated 
with hotter Hi, for which the T(Hi) = 10 4 ' 7 K photo-ionized 
model may be more relevant. This model encompasses a 
cluster of T08 absorbers with iV(Hl) « lo 134 " 13 - 6 cm -2 , 
which is near the the peak concentration of our simulated 
absorbers. Sub-Z absorbers should theoretically exist in this 
bin, although it is rare for metallicities traced by O VI to 
fall below Z. Part of this discrepancy may be explained by 
smaller 1(0 vi) (~ 50 kpc) aligned with greater i(Hl) (~ 100 
kpc), underestimating the true metallicity. COS should help 
confirm or exclude the presence of low-metallicity O VI ab- 
sorbers. 

The lowest strength aligned absorbers are slightly 
under-predicted by our simulation (1.9 Az" 1 ) relative to 
the T08 dataset (3.1 Az' 1 , a 1.1 a difference). Four of 7 
T08 weak aligned absorbers are in simple systems, and bear 
a resemblance to their proximate absorbers with iV(Hl) < 
10 130 cm -2 . The Ovi column densities are at least x3 
greater than the Hi in the data in every case, which is 
true only a quarter of the time for the simulated absorbers. 
We suggest the possibility of a non-uniform ionization back- 
ground at the Lyman limit due to the proximity of AGN 
could result in weaker Hi for an absorber, while iV(Ovi) re- 
mains unchanged, because the 8.4 Rydberg ionization field 
is more uniform as these photons have a much longer mean 
free path. This could push a number of aligned absorbers 
from the middle H I bin into the lowest bin in Figure 1111 

We find 5.1 absorbers per Az where Af(Hl) < 
10 13 4 cm -2 and £W(Ovi) = 10-30 mA, whereas there 
are none observed by T08. COS should be able to access 
such absorbers, which should indicate the enrichment level 
of the IGM at overdensities of 5-10. 

Considering strong O VI absorbers (N(0 vi) ^ 
10 14 ' 5 cm" 2 ) with aligned Hi, we find 14 (0.4 Az' 1 ) while 
T08 finds none. These absorbers are predominantly photo- 
ionized, but excluded by our ionization models, because the 
metallicity would have to be at least 10Z in most cases, 
which rarely occurs. Most of these absorbers fall off the 
7V(H i) — nH relation from Equation [12] indicating they they 
are part of a complex system with significant peculiar ve- 
locities, and possibly chance alignments. Finding no such 
absorbers in the data, we may be overestimating the Ovi 
columns of such strong absorbers as we discuss in t|5.3l 

To summarize, an in-depth examination on the 7V(H i) — 
N(0 vi) plane shows we can explain some trends in the data 
with our photo-ionized multi-phase explanation, most im- 
portantly that O vi appears t o rise quite gradually w ith in- 
creasing Hi (cf. Figure 1 of iDanforth fc Shulll l200Sh . Our 
absorbers however appear more clustered at iV(Hl) = 
1Q i3. 5-14.5 cm -2 compare d with the T08 dataset, and we 



believe a varying ionization background near 1 Rydberg is 
the best way to distribute over a wider range of H I column 
density. We encourage the use of Hubble/ 'COS to greatly ex- 
pand the aligned H I — O vi weak absorber statistics, as these 
are one of the most straight-forward and effective ways to 
constrain ionization conditions of metals and the enrichment 
levels of filamentary structures (5 ~ 10). 

4.4 Turbulence in the IGM 

By forwarding the d32n256vzwl50-bturb model, we are 
claiming that there is some form of turbulence in the low- 
z IGM well below the resolution of our simulation that is 
responsible for broad O vi lines. We arrive at this claim by 
exploring every other broadening mechanism (spatial, tem- 
perature, and instrumental), and finally settling on sub-SPH 
motions as the only viable explanation. We use "turbulence" 
as a blanket term to cover velocity shear, bulk motions, 
shock disturbances, and other random motions unresolvable 
in our simulations. Our main point is that such motions in- 
creasingly dominate the line profiles for stronger low-z O VI 
absorbers. In this section we connect (circumstantially, at 
least) turbulence to the outflows that enrich and stir the 
IGM. We begin by showing that examples of turbulence 
in the IGM and halo gas are common , and then use sim- 
ple physical relations from iKolmogorovl (|l94li ) turbulence to 
broadly motivate our prescription of sub-SPH turbulence. 

Turbulence or some other non-thermal broadening is 
com monly invoked to exp lain O vi line widths in the IGM 
(e.g. IDanforth et al1l2006h . Aligned Hi — Ovi absorbers al- 
low a constraint on how much of a b-parameter is ther - 
mal and non-thermal. Both T08 and lThom fc Chenl l|2008bh 
agree that such systems almost always have T < 10 5 K, with 
non-thermal broadening significant if not dominant. Even in 
our models without turbulence, the same exercise reveals a 
significant non-therm al broadening, a findin g also confirmed 
in the simulations of iRichter et al.l l|2006l ). However, such 
alignments do not necessarily imply turbulent motions, and 
instead may only indicate that metal absorbers arise from 
different gas than Hi. 

A more relevant measurement may be aligned H I — 
Hell components, which do not re ly on the m e tal di stri- 
bution, and have been shown by IZheng et al.l (|2004h to 
indica te the IGM is dominated by turbulence at z = 
2 — 3. iFechner fc Reimers! l|2007h downplay the amount of 
turbulence in the IGM, finding that only 45% of such 
aligned components favor turbulent broadening. Our in- 
vocation of turbulence requires the association with out- 
flows and metals; therefore we disfavor turbulence in the 
low-densit y, diffuse IGM , which is likely unenriched. Fur- 
thermore, iDave fc Tripd (120011 ) show that the low-z Lya 
forest is well described by thermal line widths alone, and 
we reaffirm this point by showing in i|3.1l that lines with 
AT (Hi) < 10 14 cm -2 are unaffected by our addition of tur- 
bulence. 

A more direct way to measure turbulence is to look 
for velocity d ifferen ces on the smallest scales possible, as 
iRauch et al. I <|200ll ) did in lensed quasars. By observing 
adjacent C IV profiles in paired lines of sights, they find 
Sv ~ 5 kms -1 on scales of 300 pc at (2) ~ 2.7. They 
apply the Kolmogorov steady-state assumption whereby ki- 
netic energy, ^(vJce), injected at a rate eo cascades through 
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increasingly small eddies at e ~ eo until viscosity transforms 
this energy to heat. The dimensionless nature of Kolmogorov 
turbulence allows the application of a simple scaling between 
size, I, and rms velo city, f rms , su c h tha t t; 2 ms ~ (el) 2 / 3 . 
Therefore, similar to iRauch et all (|200ll ). we use the en- 
ergy transfer rate e « v 3 ms /l to parametrize turbulence. For 
Ovi absorbers we assume « rms ~ b turb and the transverse 
I ~ 1(0 vi), such that 



12.5 13 13.5 14 14.5 15 



u turb 

1(0 vi) 



(16) 



This application of transverse Sv's to 6-parameters makes 
more sense for O vi, because O VI has been shown to be dis- 
tributed on scales much larger than a kpc whereas lower 
ionization species (e.g. Civ, C hi . Si iv) show sub-k pc struc- 
ture indicating small cloudlets (|Lopez et al.ll2007t ). There- 
fore, a single broad O vi profile may comprise many small 
metal-enriched cloudl ets traced by l ow io nizat ion species of 
the so rt suggested bv lSimcoe et alj (|2006l ) and lSchave et al.l 
(|2007l 1. In this scenario O vi traces more extended metals 
at lower densities where turbulence on scales of 50-100 kpc 
dominates the line width, while small, dense cloudlets traced 
by low ionization species, possibly in sub-kpc clouds, show 
less broadening due to the Kolmogorov scaling. The result is 
the type of low-z system commonly observed, where low ion- 
ization species have multiple thinner profiles along w ith one 
large broad Ovi profile (e.g. iThom fc Chenl l2008bh . which 
we suggest may be primarily photo-ionized. This is only a 
hypothesis, and should be tested in galactic-resolution simu- 
lations where turbulence and the formation of dense clumps 
can be better resolved. 



IRauch et ail (|200lh calculates 



10" 



for 



high-2 C IV absorbers, and uses the dissipation timescale, 



Tdi. 



< "KE > 

2e 



(17) 



to find that the turbulence should dissipate in ~ 10 yrs. The 
average velocity of the kinetic energy input, (uke), is specu- 
lated to arise from some sort of galactic-scale feedback. This 
short timescale, along with C IV arising from overdensities 
enriched recently by SF galaxies (OD06), indicates a more 
violent environment than the overdensities O vi absorbers 
trace at low-z; therefore we consider this energy transfer 
rate an upper limit for low-z Ovi. We calculate and plot e 
as a function of N(0 Vi) in the central panel of Figure [T2] 
for both the photo-ionized and CIE cases by assuming bturb 
from the top panel and 1(0 vi) from the bottom panel of 
Figure [7] 

If we use photo-ionized absorber lengths below 
10 14 ' 5 cm -2 (I (Ovi) ~ 50 - 100 kpc), we find energy trans- 
fer rates at most ~ 2 x 10 -4 cm 2 s -3 dropping rapidly at 
lower iV(Ovi). /(Ovi) grows rapidly above 10 14 ' 5 cm -2 for 
photo-ionized absorbers reducing e, but the absorber length 
from O VI in CIE compresses such absorbers to be <C 100 
kpc, so that these absorbers (often tracing halo gas) are 
able t o fit within a halo, e grows to match the lRauch et al.l 
(|200ll ) value for the strongest O vi absorbers. Such ab- 
sorbers we will argue in EI5.3I may be analogous to the 
O vi absorbers associated with the M W halo and thick disk , 
which must be collisionally ionized ( Sembach et al.l 120031 ; 
iFox et aI1l2004l ). Considering the thick disk absorbers with 
60 kms -1 corresponding to a structure on at most several 
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Figure 12. The 6-parameters from the d32n256vzwl50-bturb 
model are plotted as data points with 1 a dispersions. b b s is 
a linear fit to the T08 6(Ovi) — AT (Ovi) relation, and b noturb 
is the fit to the d32n256vzwl50 model. From these fits, the re- 
lation for b tur i, is calculated in Equation [5] If we consider the 
photo-ionized and collisionally ionized absorber sizes in the bot- 
tom panel of Figure [7] we can calculate the energy dissipation 
rate e (ce nter panel), w hich indicates lower rates than observed at 
high-z by I Rauch et al.l 11200 If ) for photo-ionized absorbers tracing 
the diffuse IGM. The median ages are plotted (bottom panel) to 
show that more turbulence may be associated with metals ejected 
more recently. 



kpc (|Savage et al.l 120 03:) , we subtract off the temperature 
broadening (17.7 kms" 1 at T = 10 5 ' 45 K), the instrumen- 
tal broadening (~ 20 kms -1 for FUSE), and safely assume 
no spatial broadening to arrive at b tU rb ~ 54 kms -1 . If we 
assume pure Kolmogorov turbulence and i(Ovi) = 3 kpc, 
we find a much higher transfer rate, 1.7 x 10~ 



Similar size assumptions for the ISembach et al.l l|2003h ab- 
sorbers, which are more difficult to estimate due to their 
less constrained distance, give b tU rb ~ 30 kms -1 and e = 
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3 x 10 3 cm 2 s 3 . These values are higher than in the IGM at 
z > 2 and with in MW molecular clouds, e ~ 4x 10~ 4 cm 2 s~ 3 
(|Larsonl[l98ll ). although more similar to H n regions e ~ 
2.5 x 10" 2 cm 2 s -3 for a 100 pc-sized region (|0'Delllll99lh . 
The latter do not show a Kolmogorov spectrum, but may be 
more related to outflows driven by supernovae, which may 
include IVCs and HVCs in the MW halo. The key point here 
is that the turbulence we invoke for most O VI systems is a 
small fraction of the energy transfer rate seen at high-z and 
smaller still compared to turbulence associated with active 
regions within our Galaxy. 

Do our absorbers follow a Kolmogorov spectrum? 
We should see evidence of an increasing b no t U rb in the 
higher resolution dl6n256vzwl50 simulation compared to 
d32n256vzwl50 if it is. Our added turbulence, bturb from 
Equation [3] should be reduced by a factor of i 1 ^ 3 ~ 0.8 
when halving the resolved length scales according to Equa- 
tion [16] if e is constant; however we see no evidence of an 
increasing b notU rb with resolution. Therefore, either turbu- 
lence begins at smaller scales with a larger e, or it does not 
follow a Kolmogorov spectrum. There may be a precedent 
for the former as simu lations of SNe-driven ISM turbulence 
bv lJoung et al.l (|2008T ) find that energy injection occurs over 
a broad range of scales. The turbulent energy injection of 
outflows may very well occur on scales less than the resolu- 
tion of our simulations (~ 20 — 50 kpc for halo gas) meaning 
that e is larger. 

Why would turbulence decline toward lower overden- 
sity? As a possible answer, we consider the ages of indi- 
vidual SPH particles defined by the last time they were 
launched in a wind. We define the average "age" of an ab- 
sorber as an O VI- weighted age of its contributing SPH parti- 
cles ( equation 0. Plotting the median age of O VI absorbers 
with la dispersions (bottom panel of Figure [T^J) shows that 
younger absorbers have higher column densities, which in 
turn have more turbulent broadening in our model. The av- 
erage age of metals in a 10 13 cm -2 is about 6 Gyr while for 
a 10 14 ' 5 cm -2 , the age is 2.5 Gyr. Note that this is an average 
age, and higher column density absorbers may be composed 
of a number of SPH particles with a range of ages, a sig- 
nificant fraction of which are below 1 Gyr. Our turbulence 
model is therefore consistent with the idea that turbulence 
is injected by outflows, and then dissipated on a timescale 
comparable to (or perhaps somewhat shorter than) a Hubble 
time. 

Finally, we stress the need for high-resolution simu- 
lations to trace outflows as they move from the galactic 
ISM into the IGM. The largest contributor of ISM tur- 
bulence in star-forming ga laxies appears to be supernovae 
l|Mac Low fc Klessenir2004 ). with > 90% of the kinetic en- 
ergy from^Ne^rivOT^urbjilence contained shortward of 200 
pc l|Joung k, Mac LowH2006r ). Although cosmological simu- 
lations are far from capable of resolving such scales, an im- 
plementation of a sub-grid model of turbulent pressure ap- 
pears conc eivable dJoung et al.ll2008l ). Meanwhile, the sim- 
ulations of iFuiita et al.l ( 20081 ) track the formation of shell 
fragments by Rayleigh- Taylor instability thought to be re- 
sponsible for Nal absorpti on tracing wind feedback around 
star- forming galaxies (e.g. iMartinl l2005aT ). The differential 
velocities produce Nal line widths of 320±120 kms -1 , which 
we would describe ascribe as turbulent broadening using our 
broad definition. We hypothesize that the hot medium driv- 



ing the shell fragments that escape into the IGM evolve into 
the low density metals traced by O VI with denser cloudlets 
traced by low ionization species (e.g. Cn, Cm, & Sill). Ex- 
tending galactic-scale simulations, including treatment for 
small-scale turbulence (< 200 pc), to follow winds into the 
IGM under an ionizing background should help show how 
turbulence affects the line profiles of various species. 

To summarize, "painting on" turbulence post- 
simulation should rightfully be considered controversial. 
However, it also likely to be unrealistic that gas on 20-100 
kpc scales pumped by galactic outflows has no internal 
motions on smaller scales. The energy dissipation rate of 
the required turbulence at low- z is just a fraction of that 
observed for high-z Civ, which traces more recent SF ac- 
tivity. Our discussion here emphasizes that IGM metal-line 
absorption profiles can be dominated by small-scale veloci- 
ties, analogous to the line profiles of Galactic H n regions, 
molecular clouds, IVCs, and HVCs. In fact, the diffuse 
Lya forest may be a rare instance of something relatively 
unaffected by turbulent velocities. Resolution appears to 
resolve better some of the small scale velocities in the very 
strongest absorbers, but if the dl6n256vzwl50 simulation is 
any guide, the resolution required to accurately model such 
small-scale motions is likely much higher. Lastly, we have 
added turbulence as a function of density, which is only a 
first-order approximation for the more relevant properties 
of environment and ages of metals; an interesting future 
study (particularly in comparison to COS data) would be to 
see how exactly turbulence originates and what parameters 
best describe it. 



5 ORIGIN AND ENVIRONMENT OF Ovi 

The primarily photo-ionized nature for Ovi absorbers sug- 
gests this ion traces metals in the warm diffuse IGM. We 
consider here the origin of various absorbers by determining 
when their metals were last injected by winds. The age- 
density anti-correlation alluded to in the previous section 
naturally segues into a discussion of environment, where we 
consider an absorber's relation to galaxies responsible for en- 
riching the IGM. We then discuss collisionally ionized O VI, 
noting that it appears primarily within halos of ~ M* galax- 
ies which we suggest is analogous to O VI associated with the 
MW galactic halo. Finally, we show simulated COS obser- 
vations of O VI absorption systems, demonstrating how such 
observations can provide insights about environment and the 
evolutionary state of IGM gas. 

Figure [13] illustrates various physical and observational 
properties of the IGM around four selected galaxies ranging 
over two orders of magnitude in stellar mass. These figures 
show zoomed-in views of some of the trends seen in Figure 
[T] Metals follow galaxies, and O VI in particular traces fila- 
mentary structures. Despite a large range in stellar mass in 
the left three columns, the extent of O VI is not qualitatively 
different. Furthermore, the stronger O VI in the group envi- 
ronment on the right panels are found around galaxies on the 
outskirts of this group and not around the M, = 10 11 ' 7 Mq 
in the center. Much of this has to do with the majority 
of O VI in our simulations tracing the diffuse, photo-ionized 
IGM rather than hot halo gas. In this section we will con- 
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Figure 13. The IGM environments in 8 X 8/i — 1 Mpc X 25 kms -1 slices around 4 different-sized galaxies at z = 0.25. Top row: Galaxy 
locations within the slice are indicated by colored dots corresponding to the galaxy stellar mass; the greyscale indicates gas overdensity. 
Second row: Colors indicate the enrichment level of the gas with the greyscale. Third row: H I column densities. Fourth row: O VI column 
densities. Bottom row: Average gas temperature of the IGM. Four different galaxy environments are displayed ranging from sub-L* (left 
2 columns), an L* galaxy analogous to the Local Group (third column), and a massive group environment (last column). The strength 
and extent of Ovi does not necessarily correlate with the galaxy mass, unlike Hi. Ovi is mostly photo-ionized in these simulations 
tracing overdensities of 10-200. Often photo-ionized Ovi is co-spatial with the WHIM and hot halo gas (regions of red and yellow). 



centrate on the locations of galaxies (top row) compared to 
Ovi absorbers (4th row). 

To quantify trends in origin and environment, we have 
modified our specexbin spectral generation code to assign 
the stellar mass (M,) and distance {r ga i) to the galaxy 
identified via Spline Kernel Interpolative Denmajjfl with the 
greatest dynamical influence for each SPH particle. Note 
that the d32n256vzwl50 simulation resolves galaxies down 
to ~ 10 9 Mq. We define a neighboring galaxy to each SPH 
particle as the one with the smallest fractional virial dis- 
tance to the particle (i.e. the minimum r ga i/r V i r , where 
r V ir oc M^y"^). To calculate the dynamical mass (Md V n ) 
from the stellar mass we use the relation from lDave] (2009), 
which shows the fraction of baryons in stellar mass for a 
given halo mass as calculated from this simulation and the 



d64n256vzwl50 simulatiorjf]. specexbin tracks age as well as 
the originating galaxy mass and launch u w i n d for enriched 
outflow particles, but we limit ourselves to considering an 
absorber's age only, saving an in-depth analysis of the orig- 
inating galaxies of absorbers for future work. The neighbor- 
ing M* and r ga i we refer to here often have no relation to 
that of the particle's originating galaxy, especially consider- 
ing the processes in galaxy evolution and structural forma- 
tion have altered the environments of O VI absorbers, which 
often trace metals released into the IGM at a much earlier 
time; this contrasts to the closely related galaxy-absorber 
co nnection found in simulati ons with the same wind model 
in lOppenheimer et al.l ((2009) at z = 6 where winds appear 
to be clearly outflowing from their parent galaxies. Unlike 



6 See the relation in the left middle panel of this paper's Figure 
5 SKID; http : //www-hpec . astro .Washington. edu/tools/skid. html 1. 
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the normal version of specexbin, we only use one snapshot 
at z = 0.25 over the range of z — 0.5 — * 0, although the 
varying ionization fields and Hubble expansions are applied 
from the redshift along the line of sight; this explains why 
ages never go above 10 Gyr. 

We plot Ovi-weighted densities, temperatures, ages, 
and neighboring galaxy masses at absorber line centers in 
Figure[l4]at z < 0.5 in 30 high-quality lines of sight from our 
d32n256vzwl50-bturb model. The p — T phase space panel 
(top left) reiterates our finding that most Ovi absorbers 
are photo-ionized near T = 10 4 K with an easily notice- 
able density- iV(0 VI ) correlation. The over-plotted logarith- 
mic metallicity-weighted contours show low-z metals have 
either cooled to T ~ 10 4 K, where O VI is an ideal tracer 
of these diffuse IGM metals, or are much hotter, T > 10 6 
K, where cooling times are long. In between is the metal 
"zone of avoidance" where the rare O VI CIE absorber traces 
rapidly cooling halo gas. We provide an in-depth analysis of 
this region in <j 5 .3 1 



5.1 Ages of Ovi Absorbers 

The ages of Ovi absorbers in the bottom left panel of 
Figure [14] generally anti-correlate with density. Absorbers 
below N(0 vi) = 10 14 cm~ 2 (more often in simple sys- 
tems) usually correspond to metals in filaments often in- 
jected during the high- 2 epoch of intense cosmic star for- 
mation. Strong absorbers trace metals on the outskirts of 
halos injected at z < 1, while some of the strongest ab- 
sorbers near N(0 vi) = 10 15 cm -2 trace metals cycling in 
halo fountains within the virial radius as described by OD08 
on timescales < 2 Gyr. Histograms of ^ N(0 Vi) are drawn 
for each property with the colored histograms stacked upon 
each other corresponding to the column density range; in- 
termediate absorbers (green & orange) hold the majority of 
Ovi in the IGM. A median age of 3 Gyrs is found with a 
spread of 0.6 dex. Interestingly C IV at z — 2 shows a similar 
age spread and median (~ 1 Gyr) in ter ms of a fraction of 
a Hubble time (lOppenheimer et alj|2007h . which is not sur- 
prising considering that C IV at this redshift traces similar 
overdensities as low- 2 Ovi (e.g. Simcoe et al. 2004; OD06). 

The trend of increasing age with decreasing overden- 
sity of metals results from the outside-in pattern of IGM 
enrichment. OD08 shows that winds launched at z = 
6 — > 0.5 all travel similar distances (60-100 proper kpc), 
thereby allowing a galaxy to enrich a larger comoving vol- 
ume at early times. The increasing mass-metall icity rela- 
tionship of galaxies as th e Universe evolves (cf. lErb et all 
2006; lTremontI et alj2004h naturally results in a metallicity- 
density IGM gradient as later galaxies enrich gas at higher 
overdensities filling less comoving volume with more metal- 
rich gas. Therefore the outer regions traced by the weakest 
Ovi absorbers are tracing the earliest epochs of IGM en- 
richment while remaining relatively unaffected by later en- 
richment. O vi is by far the best UV transition able to trace 
such regions due to its high oscillator strength, oxygen being 
the most copious metal, and its highly ionized state ideal for 
tracing low densities where there are many ionizing photons 
per iorQ. Conversely, the overdense inner regions trace gas 



7 It would be fascinating to see if nucleosynthetic yields indicat- 



recycled multiply in winds often not escaping from galac- 
tic halos, as described in OD08. These regions are relatively 
rare in a volume-averaged measurement as an absorption 
line spectrum, but may harbor the strongest lines with O vi 
in CIE. 



5.2 Galaxy Environments of O vi Absorbers 

The strong correlation between density and N(0 Vi) sug- 
gests a possible correlation between M* and A(Ovi), con- 
sidering that we find more overdense environments associ- 
ated with more massive galaxies at low- 2 in our simulations 
(OD08). The right panels of Figure [T4l do show such a trend, 
but only very weakly; there is a large range of absorber 
strengths for a given M ga i- This is most easily noticed as the 
large area covered by Ovi absorbers in the M„-age space. 

Another key variable we need to consider is distance 
from a neighboring galaxy, which we plot as a multiple of the 
virial radius, r ga i/r v i r , in Figure [151 Most O vi absorbers in 
the observed range reside several virial radii from the neigh- 
boring galaxy, suggesting they are dynamically unassociated 
with this galaxy. Histograms along the side for each column 
density bin indicate weaker absorbers lie at a greater virial 
distance from the nearest galaxy, which is on average less 
massive. The stellar mass associated with O vi absorbers 
is M ga i ~ io 9 - 5 - 10 M s ~ 0.03 - 0.1 AT. These same ab- 
sorbers are among the oldest, suggesting merely inciden- 
tal association with their present day neighboring galaxy. 
The absorbers at JV(Ovi) < 10 14 cm -2 associated with 
Mgai > 10 11 Mq are similarly incidental. This is best il- 
lustrated by the chaotic structure of Ovi abso rbers in the 
right panels of Figure ITgl lGangulv" et al. I l|2008h finds a sim- 
ilar association of their wea ker absorber s (EW < 50 mA) 
with ~ Q.1L* galaxies in the lCen"fc Fand l|2006h simulation. 

Our simulation results here are worth comparing to ob- 
servations, however we note that we consider only galaxy 
stellar masses, not luminosities, and our detection limits are 
complete down to M gai = 10 9 M , or ~ 0.01M*, which 
is better than a ny survey current l y obta ins, except very lo- 
cally. Surveys bv lProchaska et"al1 (I2006T ) and lCooksev et all 
(2008) of the galaxy-absorber connection find a large vari- 
ety of environments for O vi absorbers, which could pos- 
sibly support our claim that absorbers in the observed 
range are usually not directly associated with their neigh- 
boring galaxy. The majority of our absorbers are between 
100-300 kpc from their nearest galactic neighbor, which 
compares favorably with observations finding the nearest 
projected neighbor at ~ 100 — 200 kpc away from Ov i 
jTripp et al.ll200ll ; iTumlinson et al.ll2005l : iTripp et al.ll2006l ) , 
and the media n distances of 200-2 70 kpc from 0.1 M* galax- 
ies derived by IStocke et alj (|2006T >. 

IStocke et al.l 1 20061 ) stresses the need for deeper surveys 
below 0.1L* to find the galaxies responsible for the bulk of 
the IGM enrichment at low-2. Indeed, 66% of our interme- 
diate absorbers neighbor < 0.1M* galaxies, which is even 



ing enrichment by low-metallicity or even Pop III stars are asso- 
ciated with weaker O VI lines, although this would be difficult as 
no other species has a similar transition in the UV. High-S , /A r ob- 
servations of aligned C IV using COS could possibly probe similar 
low overdensities with the trend of derived [O/C] with column 
density possibly indicating evolving nucleosynthetic yields. 
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Figure 14. Ovi absorbers from the high-quality d32n256vzwl50-bturb sample are plotted in four phase space planes considering two 
physical parameters (O VI- weighted density and temperature), the age of gas tracing Ovi absorbers (i.e. the time since the wind launch), 
and the stellar mass of the neighboring galaxy (i.e. the galaxy with the least fractional virial distance). Logarithmic contours at 0.5 dex 
steps in the p — T phase space correspond to the metallicity- weighted density (lighter contours are less). Ovi photo- ionized absorbers 
trace metals up to an overdensity (<5) of 300 where /(O vi) drops to almost zero, and O VI at higher density must be collisionally ionized. 
Few of these latter absorbers are found, because metals have a "zone of avoidance," here due to efficient metal-line cooling through this 
region. Absorbers can be followed around the four panels by eye to see how these parameters are related for various O VI strengths. A 
density-age anti-correlation exists with O VI in CIE associated with recent activity, often in close proximity with M* galaxies. Absorbers 
below iV(Ovi) = 10 14 cm -2 do not correlate generally with environment showing little dependence on the neighboring M ga i. Such 
absorbers trace photo-ionized metals released in winds often at z > 1. Histograms along the side show the summed HN(0 vi), with each 
color corresponding to the column density range indicated by the key in the upper right panel. Less than 10% of T,N(0 vi) is collisionally 
ionized. 



more impressive considering our neighbor weighting scheme 
favors association with more massive galaxies. However, the 
ages of such absorbers are almost always greater than 2 Gyr, 
often much more, indicating that such association is perhaps 
incidental. The scenario may be related to C IV absorbers at 
z ~ 2 — 3 showing the same clust ering at several comov- 
ing M pc as Lyman-break galax i es fAdelberger et alj|2003t 
120051 ), while IPorciani fc Madaul (120051 ) find that this clus- 
tering arises even if the metals C IV traces were injected at 
very high-z. The same scenario appears to be in play for 
O vi absorbers, which often have little dynamical associa- 
tion with their neighboring galaxy and instead are injected 
at a much earlier epoch. 

If O vi is not linked to its immediate enviro nment, then 
perhaps this can explain rProchaska et al.l l|2006l ) finding 6 H I 
metal-free systems with qualitatively similar environmental 



characteristics as their 6 systems with O vi along the PKS 
0405-123 sight line; the Ovi does not care about the im- 
mediate surroundings and therefore traces a wide range of 
environments. This supports the inhomogeneous nature of 
IGM metal enrichment, where Ovi absorbers trace metals 
above the mean metallicity. We also find a significant frac- 
tion of our H 1 absorbers without aligned O vi (compare third 
and fourth rows of Figure 1131) . 

As an aside, we plot the fraction of H I absorbers aligned 
with O vi as a function of column density at various 5v sepa- 
rations in Figure [TS] as long as O vi has an equivalent width 
^ 30 mA. Note that unlike O vi aligned with H I consid- 
ered earlier, only one O vi component can be aligned with a 
single Hi component. We provide this plot as a prediction 
for comparison with current and future data. We predict 
about half of 7V(H 1) = 10 140 cm -2 absorbers should have 
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Figure 15. The distance to the neighboring galaxy, shown as a 
multiple of the virial radius, is plotted against that galaxy stel- 
lar mass for all O VI absorbers. The color and size key for O VI 
column density in Figure Il4l applies here. Dashed lines show the 
physical distance assuming z = 0.25, and the colored histograms 
correspond to the number of absorbers in each O VI column den- 
sity bin. Stronger O VI absorber are more likely to be dynamically 
associated with more massive galaxies, with the strong absorbers 
almost always within 2 r„; r . The range of physical distances seem 
similar to observed impact parameters. We obtain the virial radius 
by using the relation from lNavarro. Frenk. fc White! <|l997T ) for a 
dark matter halo; we convert from stellar to dyna mical mass us ing 
the relation in the left center panel of Figure 1 of lDavel j2009l) . 



O VI within 80 km s 1 , however only 20% have well-aligned 
Ovi. At jV(Hi) = 10 15 " 16 cm" 2 , O vi lies within 80 kms" 1 



% of the time, but the simulations suggest such absorbers 
are slightly less likely to be well-aligned than the bin below 
it- a possible signature of Ovi in CIE as we discuss next. 



5.3 Collisionally Ionized O VI 

The most obvious difference between collisionally ionized 
(CI) and photo-ionized (PI) Ovi is environment. The few 
strong absorbers at CI temperatures in Figure [14] trace gas 
within or just outside halos of ~ M* galaxies undergoing 
rapid metal-line cooling in the metal "zone of avoidance." 
Their relatively young ages, ~ — 2 Gyr, indicate associa- 
tion with recent enrichment. The strongest Ovi absorber in 
PKS 1405-123 (iVfOv i) = 10 14 8 cm -2 at z = 0.1671) from 



iProchaska et ail (|2006h fits this description almost perfectly 
the Ovi appears to be CI, nearly Zq, and lies 108 kpc from 
a 4.1 L* galaxy (i.e. < r„j r ), which show s evidence for a 
burst of star formation ~ 1 Gyr ago. While IProchaska et ail 
(2006) expresses surprise no other galaxies lie near this ab- 
sorber, we find this appropriate because our similar simu- 
lated absorbers depend predominantly on the environment 
associated with the parent halo; this is a clear contrast with 
PI absorbers in the diffuse IGM. 

The strongest IGM O vi absorber yet observed at z < 
0.5, the jV(Ovi) = lO 14 ' 95 cm -2 absorber in PKS 0312- 



100 



80 



0) 

c 

tsfl 



<5v 
<5v 
<5v 
6v 
6v 



4 km/s 
8 km/s 
16 km/s 
32 km/s 
80 icm/s' 




14 15 
log[N(H I) (cm" 2 )] 



16 



Figure 16. The alignment fraction of Hi absorbers with 
Ovi with at least EW = 30 mA at various Sv's us- 
ing the d32n256vzwl50-bturb model. Bins from ./V(Hl) = 
jq12.9— 15.0 cm -2 are 3 jgj, w i(j e with the last bin spanning 
JV(Hl) = lO 15 " 16 cm" 2 . Bins above iV(Hl) = 10 144 cm" 2 
are much more likely to be aligned with O VI at some level. The 
strongest H I absorbers are associated 80% of the time with O VI, 
but are less likely to be well-aligned, which we suggest is a signa- 
ture of O VI near or in halos more often collisionally ionized. 



770, s hows multi-phase behavior according to lLehner et ail 
(2008), where a PI Lyman limit system is aligned with O VI- 
bearing gas with a likely CI origin. This work finds an 0.7 L* 
galaxy at an impact parameter of 38 kpc, with a morphology 
indica ting a possible galaxy merger. Indeed, iLehner et ail 
(2008) interpret this system as possibly representing outflow 
materials cycling in a halo fountain as the impact parameter 
is consistent with a halo origin for this system. 

We propose such CI Ovi absorbers are primarily re- 
lated to those observed in IVCs and HVCs ( e.g. lSavage et ail 
l2003t ISembach et alj|2003 ; iFox et aill2006r i associated with 
the MW halo or possibly the intragroup medium of the Local 
Group. This is consistent with the idea that CI Ovi is the 
interface between a hot, tenuous halo medium at T ^ 10 6 
and press ure-confined dense clump s of ~ 5 x 10 6 M© pro- 
posed by iMaller fc Bullockl a2004T ) , and supported by the 
observations of the ILehner et ail (|2008l ) system. The ob- 
ser ved HVC jV(Ovi) — 6( Q vi) trend is also consistent with 
the lHeckman et ail (2002) relation resulting from radiatively 
cooling flows of hot gases passing through the coronal tem- 
perature regime, yielding CI Ovi. 

One problem is that our absorbers are too strong; we 
find 14 very strong absorbers (iV(Ovi) ^ 10 15 cm -2 ) over 
Az — 35 while not one such absorber is found in any quasar 
sight lines or associated with the MW halo. Such absorbers 
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rarely occur and usually require CI O VI, because the path- 
lengths for PI O VI are too long. Our simulations are very 
under-resolved to adequately model such dense clumps us- 
ing mspH ~3x 10 7 M s . Halo O vi in CIE is likely to exist 
in many smaller structures creating weaker lines, which also 
may explain their wide b-parameters as one component is 
made up of many smaller parcels of gas at random, turbu- 
lent velocities. Hundreds or even thousands of SPH particles 
per ~ 5 x 10 6 Mq cloud are likely needed to model the large 
range of densities in pressure equilibrium. 

Another reason for excessively strong CI O vi absorbers 
is possibly overestimated metallicities, assuming they are 
similar to MW HVCs. HVCs in general show metallicities 
of 0.1-0.3 Za <|Tripp et all l2003t ISembach et ail l2004bl ; 
iFox et al. I I2005T ). with rar e exceptions of su per-solar HVCs 
near the galactic disk (e.g. IZech et"all[20 08) ; conversely, our 
halo absorbers average Zq. Our winds are confined in M* 
halos recycling in halo fountains in a cycle less than a Gyr, 
which may be too fast. Again, resolution may be an issue as 
our single SPH particles fall straight through hot galactic ha- 
los, when they physically need to be pressure-confined. Two- 
phase boundaries in SPH are difficult as a particle shares 
a similar smoothing length and therefore density with its 
neighbors while it is at a very different temperature. Over- 
estimated ram-pressure stripping reduces tangential veloci- 
ties adding to this numerical a rtifact at the resolutions we 
explore fsee lKeres et al.l ( 2008) §4.3 for a discussion of this 
"cold drizzle"). If slower recycling is more realistic, fewer in- 
stances of metal enrichment occur, leading to lower metal- 
licities while reducing our (overly high) star formation rates 
in massive galaxies. Accurate hydrodynamical modeling of 
Ovi halo absorbers requires an increase of several orders 
of magnitude in resolution as well as non-equilibrium cool- 
in g and ionization. T h e non -equilibrium ionization fractions 
of iGnat fe Sternberg (|2007T ) may alleviate the overestimate 
of strong CI O vi absorbers by allowing lower /(O Vi) at 
T < 10 K making more weak O VI halo absorbers with a 
larger filling factor. 

Two characteristics of CI O VI are their multi-phase 
nature and their mis-alignment with Hi. The multi-phase 
nature of O VI appears in observations showing coinci- 
dent C hi and strong Hi in the strongest O vi absorbers 
(Prochaska et al. | |2004 iDanforth et al. 2006:; ICooksev" et al.l 
120081 ). Often our strongest absorbers show evidence that 
O vi itself is multi-phase, since such absorbers occupy tem- 
peratures in between PI and CI Ovi (see the broad range 
of temperatures in the red histogram in the upper left panel 
of Figure Q3J). The CLOUDY tables we use require that ab- 
sorbers where 10 4 ' 5 < T < 10 5 ' 3 K and S > 100 have multi- 
phase Ovi since /(Ovi) values are minuscule in this p — T 
phase space. 

The mis-alignment characteristic appears when we con- 
sider only simulated A^Ovi) above 10 13 ' 5 cm -2 ; such ab- 
sorbers where T > 10 5 K are well-aligned (5v ^ 8 kms -1 ) 
only 36% of the time compared to 55% for lower tempera- 
ture absorbers. Very mis-aligned absorbers (8v > 80 kms -1 ) 
occur 21% of the time for CI absorbers and only 12% for the 
cooler sample. The fact that only 7% of the T08 absorbers 
are what we consider very mis-aligned, including those with- 
out any Hi at all, seems to support the case for predomi- 
nantly PI O vi. 

We note the analytical model of iFurlanetto et all (|2005l ) 



can explain all Ovi as collisionally ionized when consider- 
ing structural formation shocks propagating out to a 8x the 
virial radius. It is not straightforward to quantify how far 
our virial shocks propagate, but qualitatively from simula- 
tion animation^] it doesn't seem to propagate much beyond 
2 r v i r . Moreover, virial shocks don't carry metals as they 
propagate outwards. Finally, virial shocks at these distances 
are usually at low enough densities that the high ionization 
parameter makes O VI in CIE impossible; any O vi must be 
photo -ionized at lower temperatures. While lFurlanetto et al.l 
|2005l ) suggests an even split between CI and PI O vi if the 
shocks go to 4 r v i r , our CI O vi is rarely found outside 2 r v i r . 
Like them, we find strong absorbers should remain nearly 
unaffected by the extent of virial shocks considering none 
of our absorbers over JV(Ovi) = 10 14 ' 7 cm -2 lie outside of 

2 Tyir- 

The frequency of our O vi systems with N(0 vi) ^ 
10 14 ' 5 cm -2 (i.e. summing all components within 
100 kms -1 ) is 1.8 Ajz -1 . We find a somewhat higher fre- 
quency of intersecting halos ^ 0.1M* halos (lAAz^ 1 for 
> M* and 1.8Az _1 for 0.1 — 1M*), however the frequency 
of a strong absorber within a virial radius of such galaxies 
occurs about a tenth of the time. While we find CI Ovi 
relates to the immediate environment of a galaxy in a way 
PI Ovi does not, not all strong absorbers are in CIE and 
more often trace gas within two virial radii around sub-M* 
galaxies. A fascinating survey with COS would be to explore 
sight lines intersecting galactic halos for a range of galaxy 
masses and types (?C 150 kpc). The incidence of O VI as a 
function of impact parameter will provide a handle on the 
filling factor of O VI, which we suggest here is more often 
collisionally ionized with i n hal os. A similar survey of C IV 
at z < 1 by IChen et al.l (|200ll ) finds a sharp cutoff above 
~ 100 kpc. However, strong O VI associated with H I does 
not imply that Ovi is within the halo, as we have shown 
that the associated Ovi often arises from lower overdensi- 
ties outside. O vi may not have a similarly sharp cutoff at a 
specific impact parameter as C IV if this is the case. 

To summarize, our simulations preliminarily suggest 
that the vast majority of Ovi absorbers in quasar sight 
lines are likely photo-ionized. Our CI Ovi fraction is 
less constrained; while we find a small fraction (< 10%), 
this may increase if we could resolve the structures re- 
sponsible for CI O VI, which we feel are related to 
the HVCs in our Galaxy and necessary to explain the 
strongest O vi systems. However, if we consider that CI 
Ovi is associate d with stronger absorbers, as other simu- 
lations also find |Cen et al.ll200ll ; iFang fc BrvanlbOOll ). the 
steeper fit to the differential column density distribution 
(d 2 n/dzdN(0 vi) oc iV(Ovi)" 2 , DS08) compared to other 
species {d 2 n/dzdN{Yii) oc iV(Hi)- 1 ' 7 , d 2 n/dzdN(C 111) oc 
N(C in) -1 ' 8 , also DS08) suggests there is not a large reser- 
voir of CI O VI at high column densities. H I and C ill have 
rising ionization fractions with higher density at PI tem- 
peratures; this is a clear contrast to /(Ovi), which drops 
sharply at rising densities, possibly explaining part of the 
difference in the power law fits. Hence the CI versus PI Ovi 
argument is not completely settled, but it seems very likely 
that most weak and intermediate Ovi absorbers are photo- 
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ionized, and that collisionally ionized O VI is usually associ- 
ated with galactic halo gas. 

5.4 COS Simulated Observations and 
Environment 

As a prelude to the types of observations enabled by Hub- 
6/e/COS, we consider Ovi simulated absorption systems 
from a bright quasar continuum generated using the COS 
Spectral Simulatoijf] in Figure [T7] where we assume 10,000 
second integrations at all bands. Our point in this figure is to 
contrast the metal-line systems of PI and CI O VI absorbers 
to see how environment can lead to significant observational 
signatures. The PI Ovi absorber with 6(0 vi) = 35 kms" 1 
at z = 0.25 is well-aligned (Sv = 0.9 kms" 1 ) with a Hi com- 
ponent where AT (Hi) = 10 150 cm" 2 and 6(Hl) = 47 kms" 1 
from Ly/3. This is a simple system with well-aligned C ill 
and C IV tracing gas 195 kpc from a 10 9 ' 8 Mq galaxy (4 r V i r ) 
with an average age of 4.4 Gyr and a metallicity of 0.45 Zq. 
We consider this absorber an alogous to the z = 0.207 02 ab- 
sorber in HE 0226-4110 from iThom fc Chenl |2008bl ); their 
Hi is slightly stronger and broken into two components al- 
though it appears centered on the O VI. 

We show the density and temperature (mass-weighted) 
in the bottom two panels in physical space (blue) and veloc- 
ity space (red). Peculiar velocities dominate over the spa- 
tial distribution even in a PI absorber tracing the diffuse 
IGM. At the line center different parcels of gas overlap re- 
sulting in a multi-phase PI absorber (T(Ovi) = 10 4 ' 1 & 
T(Hi) = lO 4 ' 3 K; <5(Ovi) = 47 & 5(Hi) = 86). In this 
case, gas from surrounding voids accreting into the local fil- 
ament heats the extended gas extending ~ 1 Mpc around 
the galaxy to T ~ 10 5 ' 2 K. This gas is mostly unenriched, de- 
spite being at CIE temperatures; any oxygen present would 
be photo-ionized to a higher state at such low density any- 
way. This suggests that there may be considerable WHIM 
gas that cannot be traced by any metal lines, owing to the 
small filling factor of metals in the IGM. 

To demonstrate a typical CI O VI absorber, we show a 
10 14 ' 2 cm -2 component, which happens to be aligned with 
a weak BLA (JV(Hi) = 10 13 ' 5 cm" 2 , 6(Hl) = 50 kms" 1 ). 
In contrast to the type of simple aligned absorbers T08 ob- 
serves, which indicate photo-ionization, the Hi distribution 
here is multi-component, extended, and weaker (another 
10 13 ' 5 cm -2 Hi component lies 80 kms -1 away). This 
system is multi-phase with T(Ovi) = 10 5 ' 01 K indicating 
possibly some weak aligned PI Ovi, while T(Hl) = 10 4 ' 16 
K. BLAs usually do not primarily trace gas at tempera- 
tures above 10 s K, often because aligned T — 10 4 K gas 
dominates the absorption. Lying 270 kpc from a massive 
10 11 ' 5 Mq galaxy, this system traces gas with Z = 0.45 Zq 
in the intragroup medium just outside of the halo. Unlike 
most other CI absorbers, which are often stronger, the met- 
als traced are older, 6.6 Gyr. An observed analog may be the 
z = 0.1212 system in H1821+643, which indicates CI O vi 
|Tripp et alj|200ll ). although the associated Hi is somewhat 
stronger and broader in that system. 

As an example of extreme multi-phase behavior, we 
show another system in Figure [T8l Although the CI absorber 

9 http : //arlsrv. Colorado . edu/cgi-bin/ion-p?page=ETC . ion 



at z = 0.36 is the same column density and has nearly the 
same metallicity as the PI absorber in Figure 1171 almost 
everything else about it is different. The absorber is wider, 
6(Ovi) = 51 kms" 1 , and is paired with a component 103 
kms" 1 away (AT(Ovi) = 14.4 cm" 2 & 6(0 vi) = 47 kms" 1 ). 
The main absorber is clearly in CIE, T(Ovi) = 10 5 ' 45 K at 
5(Ovi) = 350, and is slightly more than 1 r v i r (195 kpc) 
away from a 10 11 ' 1 Mq galaxy; this indicates the absorber 
arises in the intragroup medium with ages of ~ 2 Gyr. Pecu- 
liar velocities of several hundred kms" 1 heavily confuse the 
spatial distribution, creating one of the most extreme multi- 
phase absorbers in our sample. This is a Lyman limit system, 
7V(Hl) = 10 17 ' 9 cm -2 , only because of a chance alignment 
with the ISM of a dwarf galaxy (M ga i = 10 8 ' 9 M Q ) centered 
only 5 kpc from the line of sight. In our sample of Az = 35, 
this is the Hi absorber closest to any galaxy tracing the 
highest overdensity, <5(Hl) = 5 x 10 4 , and is completely co- 
incidental. Self-shielding may raise iV(Hl) substantially, but 
almost any other line of sight passing through this region 
would intersect the more extended intragroup O vi and prob- 
ably be mis-aligned with a much weaker BLA tracing halo 
or intragroup gas. The lower ionization metal species arise 
primarily from the dwarf galaxy ISM. This absorber is simi- 
lar to the z = 0.202 8 absorber observed in PKS 0312-770 by 
iLehner eT al. (2008), although in that case the Lyman limit 
system is possibly co-spatial with the O vi-bearing gas and 
not the ISM of a separate galaxy. 

The type of simulation analysis demonstrated here for 
these three systems should provide a prelude to a future di- 
rection in using simulations for metal-line absorber model- 
ing. Hydrodynamical simulations are necessary to model the 
complex environment and peculiar velocity field, the appar- 
ently inhomogeneous distribution of metals, and eventually 
the ionization conditions resulting from spatially and tem- 
porally non-uniform ionization sources. The first challenge 
is to be able to accurately simulate the complex low- 2 sys- 
tems observed by STIS and soon COS, and the second is to 
understand the physical and evolutionary processes within 
the simulations with an emphasis on environmental context. 

Our simulations appear to be able to generate such com- 
plex absorbers, however we admit that we have to search 
over a pathlength of Az = 35 to find just the right CI 
absorbers that reproduce the properties of some of the 
most well-studied low-z Ovi systems found among lines of 
sight covering a much smaller path length. Accurate self- 
consistent hydrodynamical modeling, especially in CI ab- 
sorbers, appears to be a ways away, requiring several order 
of magnitude of improved resolution, radiative transfer, and 
non-equilibrium effects. We are more confident about our 
ability to model the PI absorbers tracing the diffuse en- 
riched IGM, but we only display one, as these are a fairly 
homogeneous population. 

As for environment, we have only taken a first step in 
that direction when considering the mass and distance to 
the neighboring galaxy and the age of the absorber. The 
characteristics of the originating galaxy, the velocity it was 
launched at, the number of times a metal has been recycled, 
and the spread in ages of metals within an absorber are just 
a few of the things we have considered here. Such consider- 
ations are more informative for photo-ionized Ovi because 
these absorbers appear to be relatively unrelated to their 
neighboring galaxy. 
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Figure 17. Two metal-line systems simulated at prospective COS resolution and S/N. Absorption lines and physical parameters 
plotted in red are in velocity space, while blue indicates the spatial distribution of density and temperature in the bottom two panels 



without peculiar velocities. The left system at AT(Ovi) 



10 1 



cm is comprised of photo-ionized gas well-aligned with H I at 



JV(H i) = 10 5 cm -2 tracing gas 195 kpc from a ~ 0.1M* galaxy. The right system is a 10 cm -2 Ovi absorber in CIE lying 270 kpc 
from a super-M* galaxy showing a weak asymmetric Hi profile and no lower-ionization metals lines. Both absorbers have metallicitics 
of 0.45 Zq. Peculiar velocities dominate over spatial distribution of metals and baryons as in the systems. The two systems, comparable 
to actual observed systems, emphasize the different environments from which O VI absorbers arise. 



6 SUMMARY 

We examine the nature and origin of O VI absorbers in the 
low redshift (z < 0.5) IGM using Gadget-2 cosmological 
hydrodynamic simulations including outflows, exploring a 
variety of wind models and input physics variations. We de- 
termine the best-fitting model by comparing to a suite of 
observed Ovi and Hi absorber statistics, including the cu- 
mulative equivalent width distribution, 6(0 vi) as a function 
of AT (Ovi), and the JV(Ovi) as a function of AT (Hi). 

Our first main result is that only a model where we 
explicitly add sub-resolution turbulence is able to match all 
observations. The observations motivating this conclusion 
are the progressively larger line widths for stronger O vi 
systems, along with the high incidence of large equivalent 
width absorbers. We discuss this further later on. 

The second main result is that the vast majority of Ovi 
absorbers are photo-ionized, with only a few strong systems 
being collisionally ionized. This result does not depend on 
turbulence or any other input physics variations. The gov- 
erning variable in photo-ionized O vi strength is density, 
which steadily increases from overdensities of ~ 10 for the 
weakest observed absorbers tracing metals in extended fila- 
ments, to ~ 200 for the strongest observed absorbers trac- 
ing metals in or around galactic halos. Increasing metallic- 



ity with density also helps O vi column densities become 
stronger. 

Our third main result is that metals in the IGM are 
distributed inhomogeneously, and that Ovi absorbers prefer- 
entially arise from over-enriched regions. The average ab- 
sorber traces metals 5x above the mean metallicity-density 
relationship in our simulation (0.15 — TO Z©). Only 1.3% 
of the IGM volume enriched to greater than 0.1 Zq, with 
a filling factor of only 0.3% resulting from weaker winds 
working nearly as well. We also test a uniform distribu- 
tion of IGM metals, finding that this creates too many weak 
absorbers. Our results imply that applying the mean IGM 
metallicity-density gradient to gas everywhere in order to 
study metal absorption is not appropriate; the scatter in the 
global metallicity-density relationship critically influences 
Ovi absorbers. 

The clumpy metal distribution is an important reason 
why most Ovi is at photo-ionization temperatures, because 
it strongly enhances metal-line cooling. Oxygen is a power- 
ful coolant and is (not coincidentally) particularly efficient 
in the narrow temperature regime where Ovi has a colli- 
sional ionization maximum. This creates a "zone of avoid- 
ance" for enriched IGM gas between 10 4 ' 5 - 10 5 5 K. Previ- 
ous theoretical studies that did not account for metal line 
cooling incorrectly predict more collisionally ionized sys- 



30 B. D. Oppenheimer & R. Dave 



Collisionally Ionized VI (10 145 cm- 2 ) 




i i i i i i i i t i i i i 

200 200 

v (km s _1 ) 

Figure 18. Centered on a 10 14 ' 5 cm" 2 Ovi component in CIE, 
this multi-component system is a result of O VI in the intragroup 
medium 195 kpc from a M* galaxy aligned with a Lyman limit 
system AT (Hi) = 10 17,9 cm -2 arising from a chance alignment 
with the ISM of a nearby dwarf galaxy, M ga \ = 10 8 9 Mq, which 
is also responsible for the strong C III and C IV. Overlapping gas 
in velocity space aligns Hi gas at 10 4 K with Ovi gas at 10 5 K. 



tems. Non-equilibrium ionization allows collisionally ionized 
O VI to extend to lower temperatures l|Gnat fc Sternberg! 
l2007h . but the differences are only important for relatively 
rare high density regions, and do little to alter the over- 
all Ovi statistics. The hig h fra ction of aligned ab sorbers in 
the lThom fc Chenl (|2008af ) and lTripp et all (|200St ) datasets 
supports most O VI absorbers being photo-ionized. While ob- 
served complex, multi-phase systems with mis-aligned com- 
ponents indicating CIE Ovi may be under-represented in 
our simulations, it is unlikely that would dominate typical 
Ovi systems under any circumstance. A key implication of 
our findings is that Ovi is generally a poor tracer of the 
Warm-Hot Intergalactic Medium (WHIM), and cannot be 
used to infer the WHIM baryonic content. 

In depth study of the properties of O VI absorbers 
and their associated Hi reveals that multi-phase temper- 
atures and densities are often needed to explain their 
properties. Metal-line cooling allows inhomogeneously dis- 
tributed O vi to cool slightly below H I- weighted tempera- 
tures, while stronger Ovi absorbers often trace lower over- 
densities than aligned strong H I, which typically traces gas 
inside halos. Such a multi-density photo-ionized model ap- 
pears necessary to explain strong O vi absorbers aligned 
with JV(Hi) > 10 15 cm" 2 w hen considering ab sorbers on the 
JV(Hl)-jV(Ovi) plane. The lTripp et all (|200ot ) dataset finds 



more weak H I absorbers aligned with O vi than our simu- 
lations, possibly indicating local fluctuations in the meta- 
galactic photo-ionizing background around 1 Rydberg. 

We find good agreement with observations for the dis- 
tribution of velocity separations between Ovi and Hi ab- 
sorbers. The typical length of a photo-ionized O vi absorber 
is 50 — 100 kpc, although the peculiar velocity field often 
dominates resulting in mis-aligned H I components, espe- 
cially those tracing higher densities. Only 42% of our strong 
absorbers N(0 vi) ^ 10 14 cm" 2 are aligned within 4 kms" 1 
while 91% are aligned within 80 kms" 1 . This indicates a 
clumpy distribution of O vi that does not exactly trace the 
smoothly varying gas in the Lycv forest, but still arises in 
the same underlying large-scale structure. 

A clear anti-correlation exists between the density 
traced by O vi and its age as determined by the last time the 
metals were launched in a wind. This trend arises in our sim- 
ulations because of the outside-in pattern of metal enrich- 
ment; outer, lower overdensities are enriched first when the 
Universe is young and physical distances are small, and in- 
ner, higher overdensities are progressively more enriched by 
later winds extending a smalle r comoving distance. This is a 
consequence of our finding in lOppenheimer & Dave (2008) 
that winds travel similar physical distances relatively inde- 
pendent of redshift or galaxy mass. While a majority of met- 
als falls back into galaxies, our work here shows that low-z 
Ovi absorbers below JV(Ovi) ~ 10 14 cm" 2 provide one of 
the best ways to observe the oldest metals that remain in 
the IGM. We suggest the fascinating possibility of deriving 
nucleosynthetic yields tracing the earliest stars if another 
high ionization species, such as C IV, can be observed at low 
densities with the increased throughput of COS or future 
facilities. 

Finally, we study the galaxy environment around Ovi 
absorbers. Photo-ionized O vi usually has little to do with ei- 
ther the mass of the neighboring galaxy, averaging ~ 0.03 — 
0.1M*, or the distance, averaging > 2r V i r , because of the 
significant galactic and large-scale structural evolution oc- 
curring while these metals reside in the IGM. This helps ex- 
plain why observed photo-ionized absorbers appear to trace 
a variety of environments, with a typical distance of 100-300 
kpc to the nearest galaxy. The opposite is true of the small 
minority of collisionally ionized absorbers. These strong ab- 
sorbers are < 2r V i r from M* galaxies and have ages indicat- 
ing activity within the last 2 Gyr, and are perhaps analogous 
to Ovi observed in HVCs and IVCs in the Milky Way halo. 
The relatively young age for strong O vi signifies that these 
metals recycle back into galaxi es multiple times in what ca n 
be described as halo fountains (jOppenheimer fc Daveil2008l ). 

The most uncertain and controversial part of this pa- 
per is the claim that significant sub-resolution turbulence is 
present in O vi absorbers. We determine the amount of tur- 
bulence needed by directly fitting the observed N(0 vi) — 
b(0 vi) relation. This increases 6-parameters to observed lev- 
els by construction, while simultaneously increasing large- 
EW absorbers to the observed frequency. We attempt to 
justify this "magic bullet", as well as compare the implied 
turbulence to other instances of observed IGM turbulence. 
We point out that various authors have argued for dumpi- 
ness in IGM metals extending well b elow the mass and 
spatial resolution of our simulations l|Simcoe et al. 1 120061 ; 
ISchave et all 120071 ). If such dumpiness exists, then surely 
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those clumps must have some relative motion, so assuming 
completely static velocities below the SPH smoothing scale 
of 20 — 100 kpc is unrealistic. Our model states that strong 
absorbers can be dominated by turbulent motions, just as 
the line profiles associated with molecular clouds, Hn re- 
gions, IVCs, and HVCs associated with our Galaxy are. The 
energy dissipation required using a Kol mogorov spectrum is 
a fraction of the turbulence observed bv lRauch et al.l (|200ll ) 
in high- z C IV absorbers, which is encouraging considering 
the IGM at late times should be calmer, although it is un- 
clear whether Kolmogorov theory accurately describes IGM 
turbulence. We show that this trend is quantitatively consis- 
tent with the idea that turbulence dissipates as metals reside 
longer in the IGM. Our turbulent scenario makes the case 
that metal-line absorbers are made up of numerous cloudlets 
with lower ionization species tracing high-density clouds 
making multiple thin profiles, and Ovi tracing more of the 
low-density regions in between creating a single broader pro- 
file; such is often observed to be the structure of low- z metal 
line systems. In short, we believe that small-scale turbu- 
lence is a viable possibility, although more detailed model- 
ing is required to fully understand its implications. Apply- 
ing a volume renormalization zoom technique to cosmologi- 
cal simulations, allowing the resampling of individual haloes 
in galactic-scale simulations, provides a realistic, near-term 
possibility for following the evolution of turbulence over a 
Hubble time. 

While it may be disappointing that Ovi is a poor 
WHIM tracer and provides only a weak handle on the IGM 
metal distribution, our results open up some new and inter- 
esting possibilities for using O VI absorbers to understand 
cosmic metals and feedback. For instance, low- z Ovi may 
provide a fascinating opportunity to study some of the oldest 
IGM metals in the Universe through weak absorbers. The 
strongest absorbers appear related to the recycling of gas be- 
tween the IGM and galaxies, pr oviding a unique window into 
how galaxies get their gas (e.g. iKeres et al.i r2005. 2008]). De- 
tecting multiple ions in weaker absorbers will provide a good 
handle on physical conditions owing to the photo-ionized 
nature of O VI. Our work here is only a first step towards 
understanding the nature and origin of Ovi in the low-z 
Universe using numerical simulations, which we hope to ex- 
tend further by considering other metal species in a greater 
evolutionary context. We stress the need to model in de- 
tail the complex metal-line systems that will undoubtedly 
be observed by COS. The computational models presented, 
while state of the art, still fall well short of what is neces- 
sary, particularly for stronger collisionally ionized systems. 
We anticipate that future observational and modeling im- 
provements will shed new light on the metal distribution at 
the present cosmic epoch with all its important implications. 
This work provides a first step in that direction. 
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